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ABSTRACT 


Part II Abstract 


A tradeoff analysis has been started using computer-gen- 
erated radar simulations and mathematical geologic terrain 
models for establishing the optimum radar sensor operating 
parameters for geologic research. An initial set of mathe- 
matical geologic terrain models has been created for three 
basic landforms. Families of simulated radar images have 
been prepared from these terrain models for numerous inter- 
acting sensor, platform, and terrain variables. Selected 
images are shown to establish that computer-generated radar 
image simulation and geologic terrain modeling afford power- 
ful tools for conducting parametric optimization studies. 
Three detailed analyses using computer-generated radar image 
simulation and geologic terrain modeling are provided in the 
appendices. 


Part III Abstract 


Two research problems pertaining to synthetic aperture 
radar for geologic mapping applications have been investi- 
gated: (1) the tradeoffs between the various sensor parame- 
ters and the quantity and quality of the extractable geo- 
logic data and (2) the development of automated cechniques 
of digital SAR image analysis. Initial work on a texture 
analysis of Seasat SAR imagery is reported; the gray level 
co-occurrence matrix approach was utilized to produce re- 
sults. Synthesis of simple terrain models, imaging radar 
models, and implementation of the two as a ground truth data 
base and simulation algorithms have been accomplished. Com- 
puter-generated radar simulations are shown for combinations 
of two geologic models and three SAR angles of incidence. 


PREFACE 


A one-year pilot study has been completed by the Univer- 
sities of Arkansas and Kansas in a joint effort to define an 
optimum SAR (synthetic aperture radar) sensor and an inter- 
pretation system for global exploration or' non-renewable re- 
sources via radar remote sensing. This report relates the 
work and accomplishments of the past year at both universi- 
ties and identifies the present status of the research pro- 
gram. 

The report is divided into three parts. Part I consists 
of introductory material describing the purpose and scope of 
work performed in the pilot study, background information, 
and an overview of the program as a whole, as well as a 
statement of the actual work performed and a summary of the 
present status of the work. Part II describes the work and 
accomplishments at the University of Arkansas, and Part III 
relates those at the University of Kansas. This structure 
was adopted because the two efforts were sufficiently dispa- 
rate that an integrated approach to preparing a report was 
impractical. In addition, analysis conducted during the pi- 
lot study are reported in three appendices. 
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PART I 

INTRODUCTION 


Chapter 1 
INTRODUCTION 


The Universities of Arkansas and Kansas have completed a 
one-year pilot study to initiate the research necessary to 
define an optimum SAR (synthetic aperture radar) sensor and 
interpretation system for global exploration of non-renewa- 
ble resources. In this initial work, terrain modeling and 
radar image simulation techniques have been developed for 
use in subsequent phases of the overall program. An inti a 1 
set of textbook-example terrain models has been developed 
and simulated radar images prepared from them. Computer- 
modeling via image simulation is concluded to be a cost-ef- 
fective method for evaluating terrain variations, predicting 
results, and defining optimum sensor parameters. In addi- 
tion, a preliminary texture analysis of radar imagery was 
performed and analysis of the radar system response for geo- 
logic targets was started. 


1.1 PURPOSE 


The goals pursued in this pilot study were (1) develop- 
ment of key research techniques (i.e., terrain modeling and 
image simulation), (2) initiation of critical analyses 
(i.e., texture analysis of radar imagery and radar system 
response for geologic targets), and (3) creation of prelimi- 
nary results for simple test cases. These goals were all 
necessary to demonstrate the feasibility of using computer- 
modeling techniques to attain the dual program goals of (1) 
defining an optimum SAR sensor system and (2) developing an 
interpretation strategy for global exploration of non-renew- 
able resources using microwave remote sensing. 


1.2 SCOPE 


Because of the low level of funding for this effort, the 
research was limited to performing a pilot study in three 
primary areas of geologic modeling: (1) mineral explora-- 
tion, (2) petroleum exploration, and (3) study of igneous 
activity. The pilot study was limited to work with simple 
textbook-example cases. 


1.3 BACKGROUND 


The search for and management of non-renewable resources 
on both regional and global scales need no lengthy justifi- 
cation. The ability of radar remote sensing to contribute 
significantly to these critical efforts has been well docu- 
mented with the studies of Dellwig et al , [1], MacDon- 
ald [23, Wing [3], Wing and MacDonald [43, Dellwig et 
al. [53, MacDonald [63, and Dellwig and Moore [73. 

Geologists have dor.c -n impressive job of imagery inter- 
pretation and of development of geologic models that require 
a minimum of calibration where image recognition elements 
such as tone, texture, pattern, shadow, size, shape, and 
context are prime discriminants. These geologic models are 
much like the three-dimensional concepts achieved in the 
field, where outcrop observations provide information for 
geologic maps (i.e., the plan view) and structural cross- 
sections of the subsurface. Radar imagery enables the geol- 
ogist to analyze systemati cally the visual pattern elements 
of the terrain (i.e., drainage, vegetation, landforms, etc.) 
for geologic significance, and develop geologic models rapi- 
dly. Models ultimately provide the geologist with insight 
into the three-dimensional structure (i.e., inferences about 
Uniform type and origins, rock characteristics, and struc- 
tural complexity). 

Since 1969, commercial radars have provided imagery for 
geologic mapping in parts of South America, Central America, 
Southeast Asia, Africa, and the United States. The most im- 
pressive radar mapping program ever conducted was the Aero 
Service/Goodyear project RADAM (RAdar of the AMazon) [83. 
The entire country of Brazil, about 8.5 million sq km, was 
eventually imaged for natural resource studies. Clearly, 

radar remote sensing has proven potential in aiding global 
exploration for non-renewable resources. 

Although to date most of the proven benefits of radar re- 
mote sensing have arisen from geologic interpretation for 
non-renewable resouces, relatively little research has been 
devoted to the development of optimum specifications for the 
radar sensor, or for optimum strategies for information ex- 
traction from radar data. Surprisingly, advances in applied 
radar geologic research (i.e., the development of techniques 
that will aid in mapping of surface/tectonic features) have 
been spearheaded by the private sector. For many years only 
a limited number of radar systems were available for geolog 
ically dedicated research programs, and evaluation consisted 
of empirical observations deri/ed from analysis of imagery. 
Research was conducted with only a limited number of combi- 
nations of sensor parameters such as frequency, polariza- 
tion, depression angle, look direction, and resolution. Un- 
til now, the geometric distortions (i.e., layover. 


foreshortening, etc.) and scattering phenomena caused by the 
various geologic features (i.e., predominate landforms) 
could not be modeled, nor could their characteristic expres- 
sion in an image be realized. 

This inability to control system parameters has been 
overcome with the advent of computer modeling techniques via 
radar image simulations. Definitive studies to optimize the 
complete radar remote sensing system (including sensor, im- 
age processing, and interpretation) can now be done. Image 
simulation provides a cost-effective method for modeling of 
terrain variation, prediction of results, and definition of 
optimum sensor parameters. 

A research program is being conducted to achieve optimi- 
zation of both an imaging radar system and a total remote 
sensing system for non-renewable resource exploration. The 
research effort involves developing mathematical models of 
complex geologic structures, producing families of simulated 
radar images for these geologic models, evaluating both the 
models and simulations, and from these developing interpre- 
tation strategies and the best radar configuration for geo- 
logic exploration. This report relates the initial work 
completed in a pilot study for this research. 
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Chapter 2 
OVERVIEW 


A proposal entitled "Analysis of Geologic Terrain Models 
for Determination of Optimum SAR Sensor Configuration and 
Optimum Information Extraction of Global Non-Renewable Re- 
sources," was submitted to NASA headquarters in July 1979. 
This proposal was submitted in response to 0STA-79-A for the 
Resource Observation Division, Non-Renewable Resources Pro- 
gram. It consisted of a systematic plan for the Universi- 
ties of Arkansas and Kansas, in a joint program, to conduct 
the necesary research over a three-year period to define an 
optimum SAR sensor and interpretation system for global ex- 
ploration of non-renewable resources via radar remote sens- 
ing. This project was funded as a one-year pilot study for 
approximately 42% of the original request (NASA Contract NAG 
9-3). Despite the shortfall and its imposed reduction in 
scope for the effort, substantial progress has been made to 
date. The significant results attained in the pilot study 
are largely the result of two important factors: 

1. The originally proposed task descriptions were 
modified to establish a research effort compatible 
with reduced external fudnding. 

2. The matching contributions of the Universities of 
Arkansas and Kansas were increased significantly 
for the pilot study. 

The research performed during the past year has been to 
develop the means by which the overall program can be accom- 
plished in a subsequent three-year period. A central goal 
of the program is to determine the interpretation strategy 
and best radar sensor configuration for geologic explora- 
tion. Accomplishment of this goal requires analysis of mul- 
tiple parameter radar imagery. Past experience has shown 
that the acquisition of multiple parameter radar imagery 
over extensive control sites is not economically feasible, 
except for specific and very limited objectives. Thus, a 
major objective of the pilot study was to develop the means 
by which studies can be conducted through computer modeling 
(i.e., simulation), thereby affording the control of varia- 
bles needed to establish mathematical description of the 
terrain-sensor interaction. Terrain modeling and radar im- 
age simulations are the computer modeling tools used to per- 
form these extensive parametric and sensitivity studies. 
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The approach followed in developing the computer modeling 
capability is a logical progression through several phases 
of increasing complexity. These phases are numerically la- 
beled 1st order (i.e., deterministic models), 2nd order 
(i.e., natural terrain models), and 3rd order 
(i.e., realistic models). 

The research performed in developing the radar image sim- 
ulations and terrain modeling capabilities is specified in 
the following section, and the various phases of the com- 
puter modeling are described in a subsequent section. 


2 . 1 WORK PERFORMED 

The following tasks were performed during the pilot 
study. 

1. Develop Geologic Terrain Model s . This task con- 
sisted of identifying the landscape expression of 
key geologic features for exploration of non-re- 
newable resources, and developing mathematical 
models (i.e., deterministic or 1st order) to rep- 
resent their geometric properties such as shape, 
slope, relief, etc. 

2. Produce Simul ated Radar Images . Simulated radar 
images were produced from the 1st order geologic 
models by use of empirical and theoretical back- 
scatter functions for an extensive matrix of para- 
meters of both radar and terrain. 

3. Perform Initial Parametric Sensitivity Analysis . 

An initial family of simulated images was produced 
and analyzed for a range of parameters from each 
of the 1st order terrain models. 

The primary focus of the work performed in the pilot 
study is the development of 1st order geologic or terrain 
models. As the study progresses, complex models wil be de- 
veloped from these simple ones. They will be used to faci- 
litate extensive parametric studies. 


2.2 TERRAIN MODELS DEFINED 

Several geologic test sites have been selected where the 
surface characteristics of the terrains have been proven to 
be correlative with the subsurface structures. These sites 
represent essentially classic examples of the landform ex- 
pression of common geologic structures of perhaps several 
erosional stages. 
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Subtly expressed geologic features are considered in 
stages as 1st order, 2nd order, and 3rd order terrain mod- 
els. These terrain models are mathematical representations 
of the geometric properties such as shape, slope, relief, 
etc., of the geologic features. The terrain models are de- 
veloped in a logical progression of phases from 1st order to 
3rd order. 

1st order models are determi ni sti c . These models are 
composed of relatively simple geometric shapes that are ana- 
lytically defined. For example, a ridge, in plan view, may 
be formed in the shape of a parabola, or ellipse, to repre- 
sent a syncline or anticline. This technique can be ex- 
tended further to consider the effect on the surface expres- 
sion of different folding periods and pitch angles. The 
main features of 1st order models are that they are com- 
pletely deterministic and are used to investigate parametric 
behavior and sensitivity of specific geometric shapes. 

2nd order model s are natural terrain elevation models. 
The 2nd order model is the actual el evat ion surface for a 
region of terrain containing a well-defined example of one 
of the 1st order models (e.g., an anticline or syncline geo- 
logic structure). A typical digital elevation model 
(i.e., a 2nd order model) can be acquired from digital ele- 
vation tapes such as those produced by NCIC (National Carto- 
graphic Information Center) [93. The problem with a digital 
elevation model as acquired is the sampling frequency. The 
elevation surface is sampled once every 30 m in two orthogo- 
nal horizontal directions. This frequency is not adequate 
to support simulating SAR's having 30 m resolution. This is 
a problem which must be solved. As a first approach, fre- 
quency plane interpolation via the FFT [10] and the techni- 
que of adding zeroes to increase the apparent sampling rate 
will be investigated. A feature of the 2nd order model is 
that the entire elevation surface will be treated as though 
it were covered with a homogeneous scattering category. 
That is, each radar image simulated will portray the ground 
in the region of the feature as though it were uniformly 
covered with a single type of scattering category. This 
means that the 2nd order model adds geologic noise such as 
drainage patterns and weathering characteristics to the 
stark portrayal of the 1st order one. In essence, the 1st 
order model allows for analysis in the absence of noise, 
whereas the 2nd order model permits evaluation in the pres- 
ence of noise. 

3rd order mode l s are realistic . The 3rd order model rep- 
resents an extension of the 2nd order one to obtain a more 
faithful portrayal of the terrain. It consists of the ac- 
tual elevation surface together with the actual terrain 
cover boundaries of a region of terrain. It is, then, a 2nd 
order model with natural terrain cover boundaries added. 
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Creation of these models for specific test sites will 
involve much more effort than is needed for 1st or 2nd order 
models and will require extensive and detailed ancillary 
land use data for each site. 

The 3rd order models will be exploited simply to provide 
confidence that the simulation transformations produce valid 
replicas of what would be recorded by an actual microwave 
imaging system. They also provide a model for evaluating 
the effect portrayed in a radar image of terrain masking by 
its cover. 
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Chapter 3 

SUMMARY AND CONCLUSIONS 


Because of the low level of funding, the research was 
limited to performing a pilot study which could be done for 
simple, textbook-example cases. In spite of this limita- 
tion, important preliminary results have shown that comput- 
er-generated radar image simulation and mathematical geo- 
logic terrain modeling are valid and valuable research tools 
for obtaining an optimum SAR configuration for exploration 
of global non-renenwabl e resources. 

The detail and fidelity portrayed in a radar image depend 
upon both the sensor system parameters and the landform ge- 
ometry. The interactions between system and terrain varia- 
bles are extremely complex. Computer simulation affords a 
number of significant advantages over actual imagery for re- 
search. Chief among them is the capability of exercising 
rigorous mathematical control over both the input and the 
transformation of the input into an image. In addition, 
simulation affords a rigorous methodology for evaluating and 
comparing the results produced by different radar configura- 
tions viewing the same scene and by the same (or a differ- 
ent) radar viewing different scenes. 

A general conclusion is reached that the research started 
and work completed in the pilot study provide the groundwork 
for attaining the ambitious goal set for the overall pro- 
gram. The different approaches followed at the two univer- 
sities are intended to converge for attaining the goals of 
the program. These approaches are discussed separately. 


3.1 SUMMARY OF WORK AT THE UNIVERSITY OF ARKANSAS 


The approach followed at the University of Arkansas is a 
logical progression through several phases of increasing 
complexity. In the first phase (i.e., the work reported in 
this pilot study), 1st order geologic terrain models have 
been developed for three classes of basic geologic land- 
forms, families of radar images have been simulated from 
them for numerous different parameters, and preliminary ana- 
lyses have been conducted. These models will be extended 
and refined in succeeding phases (i.e., in subsequent years) 
from the simple, stark terrain models of the first phase to 
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more complex and realistic terrain models in the second and 
third phases. 

Preliminary analyses have shown the simulated data to 
provide valuable information for understanding the sen- 
sor/terrain interaction as portrayed visually in an image. 
Extrapolation of these analyses indicates that this approach 
using computer modeling is feasible and that the dual goals 
of the program will be attained. 


3.2 SUMMARY OF WORK AT THE UNIVERSITY OF KANSAS 

The approach followed at the University of Kansas ad- 
dresses two basic research problems with synthetic aperture 
radar: (1) the tradeoffs between the various sensor parame- 
ters and the quantity and quality of the extractable geo- 
logic data and (2) the development of automated techiques of 
digital SAR image analysis . 

The analysis portion of the work shows promise for tex- 
ture measures of SAR images using the GLC matrix approach. 
The look angle dependence of the GLC results must be inves- 
tigated because these effects must be factored into future 
use of automated texture analysis. 

The synthesis portion of the work produced geolgic fea- 
ture data bases for linear and elliptical hogback ridges. 
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PART II 


STUDIES PERFORMED AT THE UNIVERSITY OF ARKANSAS 


Chapter 4 
FOREWORD 


Real and synthetic aperture radars are used most commonly 
for remote sensor surveys. In the usual geometry, side- 
looking radars employ an antenna with its long axis along 
the velocity vector of the platform. In this configuration, 
the antenna emits a narrow beam of electromagnetic energy to 
the side. This electromagnetic energy, transmitted at a de- 
sired microwave frequency, travels at the speed of light to 
the ground. A complex interaction occurs at the ground and 
a small fraction of the energy is reradiated, or backscat- 
tered, to the antenna. The backscattered energy returning 
to the radar from a single pulse is recorded versus time 
( i . e . , amplitude versus distance) thereby forming a single 
line of an image. A series of short pulses is transmitted 
and the returning backscattered energy from each is recorded 
in sequence to build a continuous strip image pulse by 
pulse, line by line. 


The return backscattered energy is received as a function 
of the slant range distance (i.e., the propagation time) 
from the platform to the ground. It is this property which 
introduces geometric effects such as layover, range compres- 
sion, and shadow into radar imagery. Radar imagery col- 
lected at relatively small angles of incidence exhibits pro- 
nounced layover and foreshortening in high-relief areas, 
whereas imagery taken at larger incidence angles in similar 
terrain is characterized by extensive shadowing. In an im- 
age, areas of foreshortening are bright because foreslopes 
typically are compressed into a small image area. Also, the 
percentage of energy backscattered from the foreslope of a 
hill or ridge typically is greater than that from a back- 
slope. Thus, depending on the angle of incidence, fore- 
slopes are characteristically bright and backslopes are dark 
in a radar image. 


In most cases both geometric and backscatter effects are 
simultaneously present in an image. For accurate geologic 
interpretation, the desired radar imagery is that which dis- 
torts least and enhances most the landforms or terrain char- 
acteristics most indicative of geologic features. It is im- 
portant, therefore, to determine the optimum set of system 
parameters for imaging each different class of geologic fea- 
tures . 
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Propagation and backscatter are not the only effects 
important in determining the image expression of geologic 
features. They are the principal factors among many which 
contribute to the spatial presentation and visual tone of 
each feature as expressed in a radar image. Several other 
contributing factors a."e listed in Table 1. 


TABLE 1 

Principal Radar Image Factors 

1 . System parameters 

a) Frequency (wavelength) 

b) Polarization 

c) Resolution 

d) Incidence angle 

e) Fading 

2 . PI atf orm parameters 

a) Altitude 

b) Aspect angle 

3 . Terrain parameters 

a) Slope 

b) Relief 

c) Surface roughness 

d) Surface cover 


An important means of separately evaluating the contrib- 
ution of each of the interacting parameters is offered by 
computer-generated geologic terrain models and radar image 
simulations. Geologic terrain models are mathematical rep- 
resentations of characteristic landforms indicative of sub- 
surface geologic structure. Radar images are simulated from 
these models for a variety of the imaging parameters noted 
in Table 1. The simulated radar images portray the surface 
topography (i.e., landforms) for each geologic feature char- 
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acterized in a model. The simulated radar images, together 
with the geologic models, form the basis for all the subse- 
quent evaluations. Geologic terrain models and radar image 
simulations are discussed separately in the following chap- 
ters . 
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Chapter 5 

GEOLOGIC TERRAIN MODELS 


In this initial work, several 1st order geologic terrain 
models have been created. They provide, in spite of their 
unrealistic and stark appearance, ideal constructions for 
evaluating how the radar scene changes for a specific geo- 
logic feature as radar sensor parameters are changed. An 
example of such a geologic terrain (i.e., landscape) model 
and its use in radar parameter evaluation is presented in 
Appendix A C 11] . 

The geologic terrain models constructed in the pilot 
study are listed in Table 2. They con'it'tute a basic set 
for construction of more complex geologic landforms or re- 
gions . 


TABLE 2 

Basic Geologic Landforms 

1 . Basic 1 i near 1 andf orms 

a) Hogback ridges 

b) Homoclinal ridges 

2. Basic domal 1 andf orms 

3. Basic curvi 1 i near 1 andf orms 

a) Breached anticline 

b) Breached anticline/syncline 


These features are not the only ones for which 1st order 
geologic terrain models can be constructed. They are merely 
the first to be modeled in this pilot study. Each landform 
is discussed separately in the following sections. 
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5.1 BASIC LINEAR LANDFORMS 

Geologic terrain models have been created for the 
basic linear landforms Illustrated in Figure 1. These 
represent both homoclinal ridges and hogback ridges. 

Each linear landform has been embedded five times in its 
own rectangular matrix, each placement representing a dif- 
ferent look direction to the radar. This general scneme is 
illustrated in Figure 2. As can be seen from the figure, a 
given linear feature is oriented so that with a radar view- 
ing the scene from the top of the figure, the feature is 
viewed from 0, 22.5, 45, 67.5, and 90 degrees. Simulated 
radar images created from a given geologic terrain model 
thus contain the image expression of that feature for all 
look directions as captured by a specific radar. 

First order hogback ridges are treated here as stark, 
perfectly straight, symmetrical landforms which represent 
steeply dipping strata (i.e., the dip is generally greater 
than 45 degrees) exposed at the earth's surface. Where a 
sequence of steeply dipping rocks is truncated bv erosion, 
the outcrop pattern (i.e., landforms) commonly appears as 
bands or hogbacks which, on a regional basis, are roughly 
paral 1 el . 

First order homoclinal ridges are represented as asymme- 
trical landforms which represent less steeply dipping strata 
(i.e., the dip is generally less than 45 degrees) an hog- 
back ridges. There is a notable difference in tht * pness 
as well as the length of the front and back slope , 'hese 
characteristic shapes often provide the geologist », is- 
tinctive clues for determining the strike and dip of forma- 
tions , 
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Dome- shaped structures form a roughly circular to ellip- 
tical landform pattern with strata dipping away from a cen- 
tral area. rhtll structures may range in tilt f »om small 
unbreached warps a few meters In diameter to regional fea- 
tures covering hundreds or thousands of square kilometers. 
A 1st order geologic terrain model of a dome-shaped landform 
Is Illustrated In Figure 3. 


Figure 3: Basic Oomal Feature Terrain Model 
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5.3 BASIC CURVILINEAR LANDFORMS 


The first curvilinear geologic terrain model created was 
that of a simple anticline. The most complex 1st order cur- 
vilinear geologic terrain model created to date Is the anti- 
cline and syncline structure shown In Figure 4 embedded In a 
square region (i.e., a square matrix). The simple zig-zag 
patterns are characteristic landforms resulting from differ- 
ential erosion of gently folded, plunging structures. Fold- 
ing is one of the most common types of structural deforma- 
tion and these characteristic landform patterns are 
generally found in complex mountain ranges as well as less 
deformed lowlands. 


Figure 4: Ant 1 c 1 i ne/Sync 1 i ne Terrain Model 




Chapter 6 

RADAR IMAGE SIMULATION 


Radar Image simulation refers to synthesizing via a digi- 
tal computer the image which an actual radar would have 
produced if it were flown. Radar image simulation is an ex- 
tremely useful diagnostic and research tool. Simulation af- 
fords a number of significant advantages over actual imagery 
for research. Among these are the capability of exercising 
rigorous mathematical control over both the input (called a 
data base ) and the transformation of the input into an im- 
age. In addition, simulation affords a rigorous methodology 
for evaluating and comparing the results produced by differ- 
ent radars as well as the results produced by changing vari- 
ous parameters of a specific design. Simulation also allows 
separation and individual evaluation of effects which are 
inextricably coupled in actual imagery, such as backscatter 
and propagation. 

A very general approach to simulation has been developed 
at the University of Arkansas C12’J. Figure 5 is a block di- 
agram of the computer programs. The required inputs to the 
simulation programs are a data base and a backscatter file. 
As can be seen from the figure, the simulation computer pro- 
grams can produce selected intermediate and final results: 
(1) power-map image, (2) slant-range image, (3) ground-range 
image, (4) slant-range image with noise, and 
(5) ground-range image with noise. Antenna and resolution 
effects can be incorporated by processing either in the spa- 
tial domain via convolution, or in the spatial frequency do- 
main via the Fast Fourier Transform and multiplication. 
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Figure 5: Block Diagram of Simulation Programs 


6.1 BACKSCATTER 


Radar return from the ground is described by the differ- 
ential scattering cross-section, and is often called back- 
scatter. Beckscatter determines the percentage of electro- 
magnetic energy which is reridiated back to the radar from 
each scattering element within a resolution cell. A resolu- 
tion cell is considered to consist of many individual fea- 
tures or objects in relation to the scale of the radar wave- 
length, arid each of these is termed a scattering element. A 
resolution cell is defined, for a short-pulse, side-looking 
imaging radar, to be the instantaneous area on che ground 
from which the backscatter arrives at the radar antenna si- 
multaneously from each individual scattering center. 
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Variation of the intensity or shade of gray (i.e., tone) 
in radar images portrays the relative strength of the radar 
signal returned to the receiver from point to point on the 
ground. Variation of the radar signal arises from numerous 
interacting and complex causes. Two fundamental causes are 

(1) the interaction between the ground and the transmitted 
electromagnetic energy illuminating the ground (2) the geo- 
metric or propagation phenomena associated with the fact 
that radar is a ranging device (i.e., it records the re- 
turned signal versus time and thus orders features according 
to their individual distances, or ranges, to the radar). 
These causes are both discussed in a paper in Appendix 
B [133. 

Typically, the percentage of energy backscattered from 
slopes facing a radar is greater than that from level ter- 
rain and from slopes facing away from a radar. Foreslopes 
are generally depicted in an image as brighter tones than 
level ground, and brighter still than backslopes. For the 
various terrain cover types present, this difference is gov- 
erned by the backscatter trends versus the local angle of 
incidence (i.e., the angle between the local vertical and a 
line to the radar). The geometric relationships are illus- 
trated in Figure 6. 

Changes in backscatter arise from numerous sources, among 
the more important of which are (1) complex permittivity 
(both conductivity and permittivity), and (2) Roughness of 
the subsurface in relation to the radar wavelength to the 
depth where the electromagnetic wave is sufficiently attenu- 
ated. If the geometry and complex permittivity of two tar- 
gets were identical, the target having the rougher surface 
would have the stronger return. The reason is that rela- 
tively smooth surfaces tend to reflect electromagnetic en- 
ergy specularly, whereas extremely rough surfaces, called 
diffuse reflectors, tend to reradiate energy nearly uni- 
formly in all directions. Specular surfaces produce an ex- 
tremely strong return only when the angle of incidence is 
nearly normal to the surface (i.e., when the angle of re- 
flection equals the angle of incidence). Rough surfaces 
produce a relatively strong radar return in any direction. 
In contrast, if the geometry and surface roughness of two 
targets were identical, the target having the higher complex 
permittivity would have the stronger return. 

Figure 7 illustrates the general trends of backscatter 
versus angle of incidence* The four trends illustrated rep- 
resent the average backscatter response for four very dif- 
ferent types of ground: (1) specular or smooth, 

(2) slightly rough, (3) moderately rough, and (4) extremely 
rough. The surface which is smooth in relation to the radar 
wavelength produces an extremely strong return only when the 
angle of incidence is nearly normal to the surface, and the 
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Figure 6: Geometry of Imaging Radar 


strength of the return diminishes rapidly as the angle of 
Incidence is Increased. An extremely rough surface In rela- 
tion to the radar wavelength is shown to have a backscatter 
versus local angle of incidence trend which is almost a ho- 
rizontal line, consistent with the concept of producing a 
relatively strong return in any direction. In addition, the 
backscatter versus local angle of incidence trends are shown 
for two intermediate surfaces. As might be expected, the 
slightly rough surface exhibits a more specular trend 
whereas the moderately rough one exhibits a more uniform 
trend. 

The scale of roughness in terms of the wavelength of the 
incident radiation is what is important. A given surface 
may appear almost smooth at sufficiently low frequencies 
( 1 . e . , at long wavelengths) and thus exhibit nearly specular 
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Figure 7: Backscatter Trends 


reflection. As the frequency Is Increased (i.e., as 
wavelengths become shorter), the same surface will appear 
progressively rougher until It becomes extremely rough and 
thus exhibits nearly uniform backscatter [14]. The scale of 
roughness of a surface in relation to the radar frequency, 
within the limitations of this simplistic, general discus- 
sion, defines the general shape of the backscatter curve. 

The complex permittivity, for all practical purposes, 
sets the level and range of the backscatter curve. The com- 
plex permittivity is an electromagnetic description of the 
surface material and it controls the amplitude and phase of 
the reflected and refracted waves. The complex permittivity 
depends upon many physical properties of the material. A 
moderately rough surface (in relation to the radar wave- 
length) with a large percentage of water would exhibit a 
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backscatter trend like the one shown in Figure 7. A similar 
surface that is relatively dry would exhibit a similar 
trend, but it would be biased below the wet one. 

Backscatter trends in imagery can be summarized in gen- 
eral terms as follows. For a given scale of roughness and 
complex permittivity, the smaller the local angle of inci- 
dence the higher the backscatter and thus the brighter the 
image tone. 

As can be seen in Figure 5, backscatter data are used in 
the simulation programs to model the reflectivity for each 
terrain cover type present in a data base. It is these 
backscatter data which represent the interaction between the 
transmitted electromagnetic energy and each cover tvpe. 
They account for important properties both of the ground 
such as the complex permittivity and surface roughness, and 
of the radar, such as the transmission frequency and polari- 
zation. Empirical backscatter data such as those published 
by Ulaby [15] or Cosgriff et al. [ 16] are used where possi- 
ble. Theoretical backscatter models such as Fung and 
Chan [17] or Hevenor [18] are often important sources of 
data. 


6.2 DATA BASE 

The data base is a digital replica of the ground and is 
the input mechanism for getting the geologic terrain model 
into the computer simulation programs. This digital repre- 
sentation symbolizes the ground and its surface materials 
(i.e., its cover) in a grid matrix having at least four di- 
mensions, two for the horizontal location of each point, one 
for its elevation, and at least one for its cover type. 

Digital terrain data (i.e., digital representations of 
terrain elevations) are used as the spatial information in a 
data base. These data represent either actual terrain ele- 
vations or mathematical models of terrain elevations such as 
the geologic terrain models discussed in Chapter 5. Figure 
4 is a three-dimensional plot of a mathematical model used 
as a data base (it is one of the 1st order geologic models 
presented in Section 5.2) and Figure 8 is a similar plot of 
actual terrain data. 

The data plotted in Figure 4 represent the geologic ter- 
rain model of a simple anticline and syncline feature. The 
mathematical representation of the geologic feature has been 
embedded in a square matrix symbolizing a square region of 
ground. This matrix contains 512 X 512 elements. Each ele- 
ment symbolically represents a square region of 15.625 m on 
the ground (i.e., an area of 244.14 sq m). The complete ma- 
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trix, then, represents a square region of 64 sq km, each 
side being 8 km long. Each entry in the matrix specifies 
the terrain elevation at a point on the ground. The maximum 
elevation, along the crest line, was set at 100 m. In this 
way, data bases have been created from each of the geologic 
terrain models discussed in Section 5. 

The geologic terrain model is illustrated in Figure 4 in 
a three-dimensional plot. The anticline is readily identi- 
fied as the structure having the long, smooth nose and in- 
facing escarpments, whereas the syncline has an abrupt 
scarplike nose and outfacing escarpments. Slope angles have 
been set for these features from a plunge angle of 8 degrees 
(on the nose of the anticline), to approximately 60 degrees 
for the very steep escarpments, to approximately 23 degrees 
on the more gradual slopes. 

The data plotted in Figure 8 are a subset taken from a 
NCIC (National Cartographic Information Center) digital ter- 
rain tape [9] and represent the elevation surface for a re- 
gion in the mountains of Tennessee, These data were prod- 
uced by the Defense Mapping Agency Topographic Center 
( DMATC ) from the 1:250, 000-scale series of USGS (U. S. Geo- 
logical Survey) maps, specifically from NJ 16-12E. These 
data consist of an array of 512 X 512 points or pixels (pic- 
ture elements). Each pixel represents a ground spot size of 
91,44 m X 91.44 m. The scene therefore portrays a portion 
of the ground representing 2192 sq km (i.e., a square of 
46.8 km on each side) in the southern Appalachian basin re- 
gion of Tennesee. Two types of terrain are present in this 
region, maturely dissected mountains and elongated mountain 
slopes. 
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Figure 8: Plot of Actual Terrain 


6.3 SIMULATION MODEL TRANSFER FUNCTION 

For simulation purposes, a radar is modeled as the system 
illustrated in Figure 9. As can be seen from the figure, 
the transfer function can be symbolically written as 


D = Y log [ f ( h*P r ) ] + K 


( 6 . 1 ) 


where D symbolizes the density of silver grains in the film 
of an image, y is the film gamma ( i . e . , the slope of the 
D vs logE curve), ard K is a constant depending upon the de- 
velopment processing as discussed by Goodman [19]. This re- 
sult arises from convolving the system's blurring function 
(h) with the input power (P r ) 


28 


ORIGINAL PAGE Wf 
OF POOR QUALITY 


P' r = h*CF r *N 0 ) = h*P r (6.2) 


where is the return power predicted by the radar equation 
(Moore [20])* and the noise, N 0 , is a random variable having 
the appropriate probability density function. The system's 
bl urri ng fj, nction accounts for the various antenna and pulse 
length effects. 

This result is next scaled in amplitude according to the 
characteristics of the system being modeled, thereby produc- 
ing an intensity (I) as 


i = f(p' r ) (6.3) 


which is incident upon film and recorded as a density D 


D = y log(I) + K 


(6.4) 


which is identical with (6.1). 
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Figure 9: Simulation Model Transfer Function 


6.4 POWER-MAP IMAGE 

The power-map image represents, in image form, the power 
density being reradiated from each point on the ground back 
in the direction of the radar antenna. The power map has no 
actual radar analogue, but it is an important interim result 
because it provides diagnostic evidence of backscatter ef- 
fects without the perturbations of added geometric or propa- 
gation distortions. 
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6.5 SLANT-RANGE IMAGE 


The slant-range image can be considered as either an in- 
terim product or a final image. The simulated slant-range 
image is directly analogous to an actual slant-range image. 
It embodies all the geometric distortions inherent in radar. 
The energy being returned to the radar has been sorted and 
rearranged by time, or distance, in accordance with the ge- 
ometry of the radar and the location of the individual reso- 
lution cell. Geometric distortions such as foreshortening, 
layover, and shadow all are incorporated and properly dis- 
played. In essence, the slant-range image represents a map- 
ping from the spatial coordinates (x,y,z) of the ground into 
the image coordinates (r,y), range and azimuth, respec- 
tively. 


As is illustrated in Figure 5, two different versions of 
the slant-range image can be generated, with and without 
noise. In either case, the resolution, and in the noisy 
case the fading, can be processed in either the spatial do- 
main or the spatial frequency domain. Processing of resolu- 
tion and fading in the spatial domain involves more simplif- 
ying assumptions and is thus less precise than processing in 
the spatial frequency domain. As the programs are now 
structured, the spatial domain processing is faster and thus 
less expensive, and is generally used where simple statisti- 
cal models are sufficient. Processing in the spatial fre- 
quency domain is generally done where higher order statis- 
tics are required. 


6.6 GROUND-RANGE IMAGE 

The simulated ground-range image is directly analogous to 
an actual radar ground-range image. It represents a nonli- 
near mapping from the slant-range image space back into a 
ground-range (x,y) image space. All the distortions due to 
elevation (i.e., relief) mapped into the slant-range image 
are accentuated when the mapping to the ground-range image 
is done. This mapping preserves the geometric fidelity or 
orthographic arrangement of features on a level or planar 
surface such as the mean ground surface. The location of 
any point above or below the flat surface is distorted in 
the slant-range image and is further distorted in the 
ground-range image. 
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Chapter 7 
RESULTS 


Families of computer-generated radar Images have been 
prepared via the simulation progams discussed In Chapter 6 
for each of the 1st order geologic terrain models presented 
In Chapter 5. Each family of Images contains representative 
variations of the parameters In Table 1. Each Image of a 
family portrays a specific qeologlc terrain model as It ap- 
pears to a specific radar (l.e., a system representing one 
set of the parameters). 

Figure 10 Illustrates the value of using geologic terrain 
models and computer simulation. The figure contains a sam- 
ple set of simulated radar Images prepared from the breached 
ant 1 cl 1 ne/sync 1 1 ne geologic terrain model discussed in Chap- 
ter 5 (Figure 4). This sample set of images represents the 
scene Imaged by a radar having a specific operating configu- 
ration (l.e., frequency, polarization, resolution, etc.), 
from Table 1 held constant. The only variable In the figure 
Is the angle of Incidence. 

All of the simulated radar images In Figure 10 ara 
ground-range images. These images simulate the results of 
an X-band radar (i.e., one operating in the region 
5-12 GHz) being flown across the top of each image from 
left to right at an altitude of 800 km. The radar operates 
symbolically with HH polarization (l.e., horizontal transmit 
and horizontal receive) and the antenna functions over the 
terrain formed from the geologic terrain model while attain- 
ing the angle of incidence listed in the figure. 

The radar images in Figure 10 illustrate how landforms 
change appearance both with changing angle of incidence and 
with changing look direction. Presented in the figure are 
radar images simulated for 11 discrete angles of Incidence 
in the angular range from 23 to 85 degrees (angle of Inci- 
dence is illustrated in Figure 6). Note how vividly this 
presentation illustrates the changes. To a large extent the 
changes are so vivid because they are illustrated across the 
complete angular range for a single, simple, stark feature 
without the complications introduced by geologic noise 
(l.e., landscape cover variations, weathering effects, 
etc.). To a lesser degree, the impact of changing angle of 
Incidence arises because the feature contained in the scene 
is Isolated. This simple illustration demonstrates the 
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value and utility of using simulated radar images formed 
from geologic terrain models: parametric trade-off analyses 
can be made with all variables held constant (system, plat- 
form, and terrain - Table 1) except the specific ones being 
evaluated and the resultant variations are immediately evi- 
dent. 

The two detailed evaluations of radar images simulated 
from 1st order terrain models provided in Appendices A and B 
exemplify how much detailed and useful information can be 
obtained from the simulated images. Simulating radar images 
from 1st order geologic terrain models is a powerful tool 
for discerning the relationships among system, platform, and 
terrain variables as expressed visually in the product actu- 
ally employed by the end user, the radar image. 

Figure 11 illustrates a sample set of radar images simu- 
lated from the list of basic linear landforms (cf. Table 1). 
These images were produced for the same simulation scenario 
as that described for Figure 10: X-band, HH polarization, 
800 km altitude, and flying to the right across the top of 
the image. The difference in this figure is that the geo- 
logic terrain model is changed from scene to scene while the 
angle of incidence is held~at a constant 50 degrees (note 
that changing look direction is incorporated in the images 
by the way the linear landforms are oriented in relation to 
the radar). This figure, then, illustrates how a constant 
radar configuration would portray linear features as the ge- 
ometry of the feature is changed. 

Appendix C is a short paper C 2 1] illustrating the appli- 
cation of radar simulation to actual terrain instead of 1st 
order terrain models. The terrain used is that plotted in 
Figure 8. The paper shows that the simulation computer pro- 
grams produce images which are very acceptable as radar im- 
ages, and thus demonstrates the significance of using radar 
image simulation as the principal research tool in this 
study: results obtained are applicable to actual imaging 
radars. 

These sample preliminary results show that digital image 
simulation and geologic terrain modeling afford important 
advances for this kind of study. They provide cost-effec- 
tive methods for evaluating both terrain variations and ra- 
dar parameter changes, for predicting results, and for de- 
fining the optimum sensor parameters for geologic 
applications, as illustrated in Appendix A. In addition, 
they provide ways for interpreting image distortions, and 
for evaluating separate image effects such as backscatter 
versus layover, foreshortening, etc., as illustrated in Ap- 
pendix B. Further, preliminary analyses of an optimum angle 
of incidence for geologic studies indicates the best angle 
is dependent upon the terrain type being studied. As sug- 
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Figure 11: F-xamples Simulated Radar Images of Linear Laudforms 
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gested in Appendix C, for the mix of terrain types modeled 
thus far the best configuration appears to be a system hav- 
ing at least two angles of incidence, one large and one 
small. This configuration, however, is not definitive. It 
represents conclusions reached after preliminary analysis of 
a very limited number of terrain environments (i.e., cover 
and relief). 
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PART III 


STUDIES PERFORMED AT THE UNIVERSITY OF KANSAS 


Chapter 8 
PROLOGUE 


The application of imaging radar to geologic mapping has 
provided illuminating results [22]. Several successful map- 
ping programs have been conducted using aircraft sensor 
platforms, but only recently, with the flight of the SEASAT- 
A satellite, has a spacecraft sensor platform been available 
for collecting radar imagery. Though the system parameters 
of the SEASAT-A were optimized for sensing the oceans, its 
synthetic aperture radar (SAR) has provided impressive sy- 
noptic views of the earth's surface which have subsequently 
been used for geologic mapping. However, the experience 
with the SEASAT-A SAR has revealed two research problems 
whose solutions would greatly enhance the utility of future 
space missions for geolgic mapping. 

One research problem is the quantitative determination of 
the tradeoffs between various sensor parameters and the 
quantity and quality of the extractable geologic data. Ob- 
viously no one set of sensor parameters will provide optimum 
SAR images for mapping all geologic features of interest 
simply because cf their diversity. It is important there- 
fore to establish how various sensor configurations affect 
the mapping of different geologic features. From the re- 
sults of such a tradeoff analysis the SAR system design can 
be tailored to produce the best images for extracting geo- 
logic features of special interest. For example, radar pa- 
rameters such as angle of incidence and averaging of inde- 
pendent samples are closely tied to image appearance. 

The other research problem is the development of auto- 
mated techniques of image analysis. Though the synoptic 
view of the earth's surface is one of the spacecraft sensor 
platform's major advantages, it also presents significant 
problems for data reduction. That is, satellite SAR images 
which cover large surface areas now must be analyzed manu- 
ally. Recent advances in digital image processing technol- 
ogy afford the possibility of machine-aided geologic analy- 
sis of SAR imagery. Because of the diversity of the image 
manifestations of geologic features, computer analysis is 
not expected to replace human interpretation of SAR imagery; 
rather the computer should be used to ease the required 
analysis tasks. This report describes how automated texture 
analysis can be used to separate various types of macrotex- 
ture (repetitive spatial patterns, many resolution elements 
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in size) which are indicative of different geologic 
features. Once identified* pertinent statistics (e.g., pre- 
ferred directions of 1 1 nears ) would be available to the hu- 
man interpreter. 

Though the two research problems seem disjoint they are 
actually closely related. The investigation of each re- 
quires four elements for success: 

/ (X) Radar system response modeling for rough 
terrain. 

. : (2) Radar image characterization (e.g., 1st 

Analysis ■ anc | 2nd order statistics) as a function of ge- 
ologic features. 

(3) Terrain modeling for macro- and micros- 
tructures . 

(4) Implementation of radar and terrain mod- 
Synthesis els to synthesize images. Analysis of images 

to determine useful geologic information. 

Clearly, both problems require SAR modeling for rough 
terrain. Because of the many complex interactions between 
the sensing of the image and the production of geologic 
maps, the only feasible technique for establishing the de- 
sired tradeoffs is to model accurately the SAR response to 
representative types of terrain and then to evaluate quanti- 
tatively the quality of the extractable geologic data as the 
system configuration is varied. The automated data analysis 
problem is aided by the same modeling. All successful im- 
age-processing algorithms require a “world model" to de- 
scribe the statistical characteristics of the data to be 
processed; thus, for data extraction, modeling and simula- 
tion are also required. 

The second key element is the investigation of radar im- 
age properties (e.g., 1st and 2nd order statistics or tex- 
ture) as a function of geologic features or vegetation 

cover. Si gnal -to-noi s e ratio and texture are two image 
properties which are strongly linked to the types of image 
surfaces. Quantitative relationships between terrain para- 
maters and image properties can be established by either 
theoretical systems analysis or sensor simulation. If the 
relationships between image characteristics and the image's 
geologic mapping potential could be identified, it would be 
possible to establish how the sensor parameters affect the 
geologic mapping potential of SAR images. Further, if this 
cause and effect mechanism were identified, it would provide 
a basis for using measured SAR image properties to indicate 
the presence or absence of geologic features of interest. 
The establishment of the correlation between SAR image char- 
acteristics and geologic features is thus important for sys- 
tem optimization and SAR data reduction. 
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Third, to support the first two components of the re- 
search the modeling of surface features and their spatial 
arrangement (e.g., a surface feature would be a ridge and a 
spatial arrangement of ridges would form a mountain range) 
is essential. This modeling effort would provide an under- 
standing of what the radar system is attempting to sense. 
Clearly these models describing the geologic features are 
necessary for the development of automated data reduction 
algorithms because they are the phenomena which are to be 
estimated (or extracted) from the imagery. Both random and 
deterministic models are useful to describe surface fea- 
tures, Though structures of size on the order of a wave- 
length are best modeled as manifestations of random proc- 
esses, larger structures may sometimes be modeled as either 
stochastic or deterministic; both approaches are being pur- 
sued by the authors. 

This report describes the work that has been accomplished 
in the four areas outlined, specifically (1) texture charac- 
teristics (2nd order statistics) as a function of different 
geologic feature orientations, (2) radar system models for 
simulation and radar image processing, (3) simple determin- 
istic models for large-scale geologic features, and (4) gen- 
eration of images given the simple terrain models and imag- 
ing radar models. 

The discussion of results is preceded by a few comments 
on geologic feature roughness scales. The many aspects of 
the microwave-terrain interaction are not reviewed as that 
topic is treated in many reports and texts (e.g., Daily et 
al. [22], Evans [23], Reeves et al. [24], and others). 


Chapter 9 

SAR IMAGING OF ROUGH TERRAIN SURFACES 


We differentiate among three scales of geologic features: 
(1) those of size on the order of the wavelength of illumi- 
nation (lithologic), (2) those on the order of a resolution 
cell (using the SEASAT-A SAR as an example), and (3) those 
on the order of many resolution cells (landforms). Radar 
response to category 1 is discussed in a nontheoreti cal 
treatment by Daily et al. [22]. Our research is aimed at 
categories 2 and 3. For example, we base our approach to 
the automated data reduction problem on the observation that 
large-scale repetitve geologic features give rise to repeti- 
tive patterns in the associated radar imagery. Such pat- 
terns can be viewed as texture features in SAR imagery. As 
we show in Section 10.1, available texture discrimination 
algorithms can be applied to radar scenes to separate scenes 
on the basis of 2nd order statistics. 

Microwave propagation effects due to the large-scale ter- 
rain elevation variations such as shadow, foreshortening, 
layover, and local angle of incidence effects account for 
the vivid depiction of natural terrain that SAR can produce. 
These important factors have been incorporated in the ter- 
rain modeling and radar modeling. These factors are inti- 
mately linked to texture within radar imagery; we show that 
different orientations of landforms (and their consequent 
different shadow patterns) are separable by texture meas- 
ures . 


Chapter 10 

STATISTICAL CHARACTERISTICS OF SAR IMAGERY 


We have proposed to study both 1st and 2nd order statis- 
tics within radar images and to relate them to landforms, 
structural features, radar look direction, angle of inci- 
dence, and other factors. First order statistics are point 
statistics, whereas 2nd order statistics can measure neigh- 
boring point relationships. 

Though for flat terrain the speckle statistics (/i,cr) have 
been known theoretically for many years (since the late 
1800's) as interference phenomena, the point statistics have 
not been treated theoretically for particular random process 
models of terrain in conjunction with deterministic sensor 
models. We have conducted a preliminary experiment on the 
texture characteristics (2nd order statistics) of regions 
Within SEASAT-A SAR images. The results show potential for 
a texture approach to automated data reduction. First order 
statistics investigations for SAR images of rough terrain 
models represent on-going work. 


10.1 A PRELIMINARY TEXTURE ANALYSIS OF RADAR IMAGERY 

Spectral, textural, temporal and contextual features are 
four important pattern elements used in human interpretation 
of image data. Spectral features describe the average band 
to band tonal variations in a multi-band image set whereas 
textural features describe the spatial distribution of tonal 
values within a band. Contextual features contain informa- 
tion about the relative arrangement of image segments be- 
longing to different categories, and temporal features de- 
scribe changes in image attributes as a function of time. 
When small image areas within, say, a synthetic aperture ra- 
dar (SAR) image are processed i ndependent ly on a computer, 
only the tonal and textural features are available to the 
machine. 

In many of the automated procedures for processing radar 
image data from small areas, such as in crop classification 
studies, only the average tonal values are used for develop- 
ing a classification algorithm. Textural features are gen- 
erally ignored on the basis that the poor resolution of ra- 
dar imagery does not provide meaningful textural information 






for such applications. However, for many other applications 
such as the identification of large-scale geologic forma- 
tions, land use patterns, etc., the resolution is more than 
adequate to provide textural information. Indeed, in these 
applications, texture is probably the most important image 
feature. We describe a procedure for obtaining numerical 
descriptors for characterizing the textural properties of 
segments of radar images. 

The textural feature extraction algorithm we describe has 
been widely used (Rosenfeld and Troy [25], Haralick and An- 
derson [26], Sutton and Hall [27], Haralick et al . [28]) for 
analyzing a variety of photographic images. The procedure 
is based on the assumption that the texture information in 
an image block "I" is contained in the overall or "average" 
spatial relationship which the gray tones in the image "I" 
have to one another. This relationship can be characterized 
by a set of gray level co-occurrence (GLC) matrices. We de- 
scribe a procedure for computing a set of GLC matrices for a 
given image block and define a set of numerical textural de- 
scriptors (features) that can be extracted from the GLC ma- 
trices. These textural features can be used for automated 
analysis and classification of blocks of radar imagery. 

Segments of digitally correlated SEASAT-A SAR imagery 
have been processed using the textural features. The re- 
sults are presented to show that these features can be used 
to classify large-scale geologic formations. Because the 
areal character i sti cs of texture carry so much information, 
we conclude that it is important to use textural features in 
automated radar image processing schemes except in applia- 
tions where the poor resolution of the imagery does not pro- 
vide meaningful textural information. 


10.2 TEXTURAL FEATURES 

Texture is one of the important characteristics used in 
image analysis. Image texture may be viewed as a global 
pattern arising from a deterministic or random repetition of 
local subpatterns or primitives. The structure resulting 
from this repetition is very useful for discriminating among 
the contents of the image of a complex scene. A number of 
approaches have been suggested for extracting features that 
will discriminate between different textures (Rosenfeld and 
Troy [25], Haralick and Anderson [26], Sutton and Hall C 2 7 U , 
Haralick et al. [28], Galloway [29]). Of these approaches, 
textural features derived from gray level co-occurrence ma- 
trices ( GLCM) have been found to be the most useful for ana- 
lyzing the contents of a variety of imagery in remote sens- 
ing, biomedical, and other applications (Shanmugam and 
Haralick [30], Haralick and Shanmugam [31J, Kruger et 
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al, [323, Weszka et al. C33D, Connors and Harlow [34]). The 
GLCM approach to texture analysis is based on the conjecture 
that the texture information in an image is contained in the 
overall or average spatial relationship among the gray tones 
of the image. 

The 2nd order gray level co-occurrence matrix of an image 
is defined as follows. Let f(x,y) be a rectangular digital 
picture defined over the domain 
x e [0,n ), y e [Q,n h x, y e I. Let ng be the number of 
gray levels in f. ihe unnormalized, 2nd order GLC matrix is 
a square matrix F of dimension nq. The (i,j)-th entry in F, 
denoted by P^-, is a function or the image tonal values and 
a displacement vector cT = (di,d 2 ). The entries Ffj are un- 
normalized counts of how many times two neighboring resolu- 
tion cells which are spatially separated by ? occur on the 
image, one with gray tone i and the other with gray tone j. 
That is, 

Pjj = ^ I (( m ]> n -|)> |f( m i> n i) = i» ^mg*^) ” J* 

and (m 2 ,n 2 ) - (m 1 ,n 1 ) = d ( , (10.1) 

where § denotes the number of elements in the set, and the 
indices mi, m 2 and m» r\o ta * <e on integer values in the in- 
tervals (Oi, n x ),(0i, n y ). The normalized GLC matrix P with 
entries Pjj is obtained from F by dividing each entry in p 
by the total number of paired occurrences. The definition 
of 2nd order GLC matrices can be extended to include 3rd and 
higher order GLC matrices. Though higher order GLC matrices 
may be important in some applications, much of the recent 
work in texture analysis has been based on 2nd order GLC ma- 
tr i ces . 


The 2nd order GLC matrices _are computed for various val- 
ues of the displacement vector d, and features derived from 
the GLC matrices are used for classifying the contents of an 
image . 


Some of the commonly used 
the GLC matrix are: 

1) Uniformity (sum of squares) 

2) Contrast 

3) Correlation 

<r 

4) Entropy 


textural features derived from 
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(10.2a) 
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(10.2b) 
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5) Inverse difference moment 


EE<p„ ) p /|1-j| u 

1 j J 


(10. 2e) 


6) Maximum Probability : max p.. 

u ,J 


(10. 2f) 


For a variety of imagery (aerial, micrographic, and x- 
ray) the relationship among these textural features, their 
values, and what they represent in terms of visual percep- 
tion of texture are reasonably well understood. Using fea- 
tures of the form given above, Haralick and Shanmugam (Har- 
alick et al. [28], Galloway [29], Shanmugam and 
Haralick [30]) were able to classify a variety of images 
with better than 85 % classification accuracy. 


Given the past success of image analysis using textural 
features on Landsat and other imagery, it seems reasonable 
to expect similar success in the case of textural features 
for analyzing radar images. Such has not been the case un- 
til now for several reasons, mainly the limited availability 
of radar images in digital form and the use of human inter- 
pretation of radar images. With the increasing need for au- 
tomated procedures to analyze large volumes of radar data in 
digital form, textural features are expected to have an im- 
portant role in radar image analysis., In the next section 
is an example which shows that textural features can be used 
for classifying segments of radar images corresponding to 
different geologic formations. 


10.3 RESULTS 

The textural features described in the preceding section 
were used to analyze segments of SEASAT-A SAR imagery over 
Tennessee (Figure 12). Four areas within this image were 
processed using the textural features described in the pre- 
ceding section. All four areas are in the vicinity of the 
Pine Mountain thrust and are underlain by the same rock se- 
quence: Pennsylvanian shales, sandstones, and siltstones. 

However, their surface manifestations are different because 
of the nature of the stress to which they have been sub- 
jected, orientation and nature of the fractures, and ero- 
sion. Six subimages were chosen from the four areas of 
analysis. Examples of the subimages which were analyzed in 
detail using textural features are shown in Figure 13. 

The co-occurrence matrices and textural features were 
computed for each of the 24 samples. The sample image 
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Figure 12: Original SEASAT-SAR Image 






Set 1 


Set 2 


Set 3 


Set 4 


Figure 13: Typical Subimages for Four Image Categories 
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blocks were of size 200 X 200 and the Image blocks were 
equal -probabl i 1 ty quant i zed_1 nto 32 levels. Co-occurrence 
matrices were computed f or | d | >1, 2, 4, a«,<J 8, and angle of 
3 »0, 45, 90, and 135 degrees. The co-occurrence matrices 
for t^e four angles were averaged and the uniformity, con- 
trast, maximum probability, and inverse difference moment 
features were computed from the average co-occurrence ma- 
trices. A total of 16 features were used to characterize 
the textural properties of each image block. 

Scatter diagrams of the numerical values of pairs of tex- 
tural features are shown in Figure 14, and plots of average 
feature values as a function of |d| are shown in Figure 15. 
i he results displayed in Figures 14 and 15 clearly demon- 
strate the usefulness of textural features for classifying 
radar Image segments. Though the image blocks have nearly 
the same mean and variance, the textural features separate 
the image categories very well. 
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Scatter Diagram of Textural Features for Four 
Image Categories 


Sum of Square Probability 




10.4 CONCLUSIONS 


A set of features for characterizing the textural proper- 
ties of segments of radar imagery is presented. The tex- 
tural features were computed and used for identifying some 
simple geologic formations from radar imagery to demonstrate 
the potential usefulness of textural features for radar im- 
age classification. 
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Chapter IX 

RADAR SYSTEM MODELING AND RADAR IMAGE SIMULATION 


Experience of the image simulation group at the Univer- 
sity of Kansas Remote Sensing Laboratory has led to two 
basic approaches in simulation: (1) calculate a geometry- 
corrected (range perspective) map of average return powers 
(real valued), incorporate signal fading, and subsequently 
perform smoothing operations on intensity data (Holtzman et 
al. [35,36]); (2) calculate a geometry-corrected complex 
map, and perform resolution degradation and averaging of 
looks independently on the complex data. In the latter me- 
thod, additional looks can be obtained only by addition of 
images on an intensity basis; degradation of resolution on 
coherent data does not add to the number of looks (the ran- 
dom walk in space analysis is still appropriate). 

Radar system modeling for ideal (isotropically reflecting 
point) targets is the standard SAR analysis (Brown [37], 
Harger [38]). Modeling for nonideal targets and collections 
of such targets according to random elevation models has not 
been accomplished; our efforts in the on-going phase of this 
research address these aspects of the SAR-geology interac- 
tion. 

Several studies of radar image simulation have been re- 
ported (Mitchell [39], Bell [40], Beckner and Crow [41], 
Holtzman et al . [35,42,36,43]). Because of the numerous as- 
pects involved in the implied geometry problems, data base 
construction, and implementation complexities, we forego a 
lengthy discussion here. However, we have found it neces- 
sary to derive processing algorithms to operate on the 
(real) output of the simulation software for system studies 
(resolution, looks). One incoherent filter which is useful 
for treating image (intensity or amplitude) data is pres- 
ented next. Processing on the complex outputs of our simu- 
lations is based on well known techniques (Zelenka [44], 
Porcello et al . [45]). 


1 1 . 1 COHERENT AND I NCOHERENT PROCESSING FOR RADAR 

Coherent processing will be defined as any operation on 
amplitude and phase of a radar signal; incoherent process- 
ing will be defined as operations on amplitudes only or on 
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power only. Thus, azimuth focusing and range 
are coherent operations; signal processing to 
si gnal-to-noise ratio involves summing intensity 
coherent processing). 


compression 
improve the 
images (in- 


The simulations presented in this report are not gener- 
ated from a simulated signal film (which would need to be 
filtered in azimuth and range to produce a radar hologram). 
Rather, the complex radar image is calculated directly. For 
each facet in the data base (the facet size is dictated by 
the resolution of the data base) the angle of incidence, 
range, local angle of incidence and antenna weighting are 
computed and subsequently used in the radar equation 
(Moore [20]) to predict the average return power at the an- 
tenna terminals. This value of Pp is used to calculate the 
amplitude of the quadrature voltage components as shown in 
Figure 16. When this operation is completed for the entire 
data base the complex radar image is available for either 
coherent or, incoherent processing. (Implicit in the calcu- 
lation of p r are the geometry calculations to incorporate 
shadow, layover, and foreshortening.) 


The image is defined as the sum of the squares of the 
voltages out of the detector 


v D (t 0 ) = [(V,(t)*h(t)|)‘ + (V q (t)*h(t)|) 2 ] = P 0 (t 0 ) 


t=t r 


t=t r 


( 11 . 1 ) 


This value at t 0 corresponds to some image location Xo,yo. 
Only incoherent processing can be done on Vg(to). Coherent 
processing can be carried out prior to the detector. 

A fixed amount of information is available in V i ( t ) and 
Vq(t). It is well-known that from this base of information 
many different SAR images can be formed, i.e., there is more 
than one processing route and some are more suited than oth- 
ers to various earth resource applications. Given the com- 
plex radar image I(t), i.e., 


I(t) = Vj(t)*h(t) + jV Q (t)*h(t)= v 2 (t) 


( 11 . 2 ) 


one can determine the image spectrum I(w) = F[I(t)]. If 
this image spectrum is subdivided and if intensity images 
are made from the inverse transformed subdivisions, the re- 
sultant intensity images can be added on an intensity basis 
to produce a radar image of the same scene which has poorer 
resolution, but an improved si gnal-to-noise ratio (Zel- 
enka [44], Porcello et al. [45], Goodman [19]). Thus, reso- 
lution can in a sense be traded for si gnal-to-noise improve- 
ment. If only one of the images from the subdivided 
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spectrum is examined, its resolution will be poor, and no 
improvement in si gnal-to-noise ratio will have been ob- 
tained. 


If one assumes digital image processing (as is the case 
here), the image spectrum is obtained from the image by 
means of a discrete Fourier transform; the original image is 
padded to twice the original dimensions to avoid circular 
convolution effects (e.g., if the original image is 
512 by 512 samples, its Fourier transform is a 1024 by 1024 
sampling of the continuous function I (to) ) . Because I ( x ) is 
complex, the spectrum is not symetric. If that spectrum is 
divided into four spectra of size 512 by 512, then, from the 
intensity images of the inverse transforms of the spectra an 
image can be made that has N = 4 statistics (where N is the 
number of independent samples) and one fourth the resolution 
of the original image (p = p x p ). 


This technique is feasible when the desired value for N 
is small. However, it is not uncommon to specify n > 10, or 
even N > 100 for orbital imaging applications. In this sit- 
uation the simulated radar image is produced by an incoher- 
ent degradation of an N = 1 image. It can be demonstrated 
that an incoherent filter transfer function can be found to 
operate on the spatial domain, to affect the operation of 
the multiple bandpass filtering (described in the preceding 
paragraph) in the frequency domain, based on Zelenka C44]. 
In general, it is possible to degrade the resolution and in- 
crease the number of looks (assuming the entire image band- 
width is utilized) according to the relation 


where 




(11.3) 


az. 


az. 


az. 


az. 


= azimuth resolution before incoherent processing, 
= range resolution before incoherent processing, 

= azimuth resolution after incoherent processing, 

= range resolution after incoherent processing, 

= number of independent samples in azimuth before incoherent 
processing, 

= number of independent samples in range before incoherent 
processing, 

= number of independent samples in azimuth after incoherent 
processing. 
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* number of independent samples in range after 
incoherent processing. 


The image bandwidth can be subdivided more times in azi- 
muth than range (or vice versa) to improve the number of 
looks in the azimuth dimension, independent of N r An in- 
coherent transfer function to trade resolution for Qtveraging 
has been implemented and used for results reported in Chap- 
ter 12. 


11.2 DERIVATION OF AN INCOHERENT TRANSFER FUNCTION 

The purpose of this development is to derive an incoher- 
ent transfer function (Goodman [193) which performs in the 
spatial domain the continuous (scanning) mixed integration 
process described by Zelenka [443. This procedure is com- 
monly referred to as incoherent averaging. For the commonly 
needed case of large N statistics it is computationally in- 
efficient to use the spatial frequency filtering approach to 
achieve the desired point statistics. The filter described 
simultaneously degrades resolution and achieves the needed 
large N statistics in one spatial domain operation by using 
image data (no phase data), in contrast to complex filtering 
in the spectral domain. 


Zelenka [443 has shown an ideal linear system model for a 
synthetic aperture radar system as represented in Figure 17. 
Incoherent processing of the radar data to improve the sig- 
nal-to-noise ratio is commonly implemented as seen as in 
Figure 18. The complex radar reflectivity s(x) is convolved 
with the impulse response p(x) and the data are filtered by 
scanning in the frequency plane. The detection on film rep- 
resents square law detection and the finite frequency plane 
window and window velocity give rise to time averaging, the 
output of the averaging filter being the image Z(x). Zel- 
enka has derived the equivalence between the scanning fre- 
quency plane mixed integrator and an incoherent transfer 
function |h(x)r which operates on detected (intensity) 
data. The extension of this development to two dimensions 
is straightforward. The results state that 

T/2 « « 


Z(x) 


1 

T 



U(f 1 ) U(f 2 ) H p (f - B t ) * 


-T/2 -« 


j2irx(fi-fo) 

H p (f 2 - Bt) e df-j df 2 dt 

= |s(x)*p(x)| 2 *|h(x)| 2 


(11.4) 
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where U(f) ■ F fu(x)] and the bar denotes conjugation. The 
transfer function H(f) - F Ch(x) D has the form of the low 
pass equivalent of the bandpass filters H p (f). For example, 
because the bandpass filters are uniformly weighted windows, 
i.e., H p (f)» rect [f/B 3, the impulse response varies as 
s1n(x)/x for the scanning filter, and the magnitude-square 
of h(x) is a sine-squared function. Therefore, |h(x)| l§ — * . 
sine-squared response (the incoherent impulse response h. (x) = |h(x)p). 

i nc 



s(x) 



w(x) 


Figure 17: Simple Linear System Model for SAR 
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Figure 18: Scanning Mixed Integration Processing 


We consider two cases in which the original image data 
are incoherently filtered to N looks and M looks, respec- 
tively. Then from equation (11.4) one can write 

|s(x)*p(x)| 2 *|h N (x)| 2 = I N (x) (11. 5a) 


and 


|s(x)*p(x)| 2 *|h M (x)| 2 = I M (x) 


(11.5b) 
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The impulse responses h M (x) and hfj(x)have been written in 
this manner to indicate that the respective functions have 
the same analytical form, i.e., the same envelope, but the 
null-to-null widths may differ. If all the bandwidth B of 
the complex radar image spectrum is employed. 


B - MB m - NB n 

Let g N (x) = |h N (x)j 2 and g M (x) = |h M (x)| 2 . 
(11.4) and defining I(x) = |s(x)*p(x)| 2 , 

X M (x) = I(x)*g M (x) 

I N (x) = I(x)*g N (x) 


( 11 . 6 ) 


. Using equation 
one can write 


(11.7) 


We note that because Im(x), In(x) represent image data they 
are strictly positive functions, as must be gjyj ( x ) and g^x) 
(they are sine-squared impulse responses). Next we consider 
the filtering operation as shown in Figure 19, in which we 
seek to generate Im(x) from I ( x ) using a filter with im- 
pulse response g(x), which also must be positive everywhere. 
Because the bandwidth of the data Ifl(x) is greater than that 
of I M (x) (g(x) is a lowpass filter), g(x) may exist for any 
situation in which M > N. 

Thus, we seek g(x) to solve the problem 


I(x)*g N (x)*g(x) = I M (x) * I(x)*g H (x) (11.8) 

Taking the Fourier transform of both sides of equa- 
tion (11.8), we obtain 


I(f)*G N (f).G(f) = 1(f). G M (f) (11.9a) 

or, equivalently, 

G M (f) 

G ( f ) = (11.9b) 

Because g^x) and 9m(x) are both sine-squared impulse re- 
sponses, Gf*|(f) and G|\j(f) are both triangle functions (and 
zero phase functions) as shown in Figure 20 and can be de- 
scribed as 
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(a) 



(b) 


Figure 19: Inverse Filter-M Looks from N Look Image Data 

V f > - A m (1 - |f|/B M ) |f| < B m 


* 0 

6 N (f) * A n {1 - |fl/B N ) 


elsewhere 

Kl < b n 


= 0 


elsewhere 


Therefore, 6(f) 

\ 
7C 


8(f) - 3* 


has the form 

i - m/B„i 

^TF|7B N _ 


given by 

m < b m 


elsewhere 


(11.10) 


59 





ommr&. puzs w 

OF POOR QUALITY 


Figure 20: The Triangle Spectra of G M (f) an G M (f). 


Because Bm > B^ and f is always less than over the 
nonzero region of G(f), the denominator of G(f) can be ex- 
panded in a convergent geometric series 

G(f) = 0 " (’ + ^ + ^ + ■ ■ ■ ) 

= /, . Mini (i + mn + 
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Expanding 6(f) and truncating the product beyond the 
quadratic terms yields 


G( f) iHudi |fi e 


(nag) 


The truncation introduces negligible error, for example, 
when N = 1 and MHO. When N increases, so must M to keep the 
approximation valid. G(f) can be rewritten as 


6(f) 





K I < b m 


(naa) 


where 


K t - r 


v 




Then it is fairly straightforward to find the Fourier trans- 
forms of the triangle portion of the spectrum and the para- 
bolic portion, with the aid of the scaling rule (Brace- 
well { 10!1 ). The impulse response of the approximation in 
equation (11.13) was derived as. g(x). 
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where 


2 n 


4it fc Bx 

| <i nr**u . n»r 

V2(MN-r) 


The original requirement on g(x) was that it he positive ev- 
erywhere, so application of the impulse response must be 
preceded by verifying that g(x) is non-negative. This was 
done experimentally for many values of M and N to examine 
the series (equation 11.11) truncation effects. 
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11.3 APPLICATION OF THE INCOHERENT FILTER 


The importance of a technique for Incoherently degrading 
an amplitude or intensity radar image (in tt, spatial do- 
main) to an arbitrary number of looks stems from the fact 
that few investigators have access to the complex radar 
data. Likewise, it is not always feasible to generate the 
complex radar simulations described before, which can be 
filtered in the frequency domain to increase the number of 
looks at the cost of resolution. 


Chapter 12 

RADAR IMAGE SIMULATION: PRELIMINARY RESULTS 


The purpose of this section is to present a sequence of 
radar image simulations which illustrate the simulation 
model and its successful implementation. Also these results 
show that simulation can be used to investigate the complex 
interactions (tradeoffs) between the terrain and sensor's 
images. In this preliminary tradeoff study, seven terrain 
features were modeled: linear hogback ridges oriented at 
five different angles with respect to the radar look direc- 
tion and two symmetrical hogback ridges following elliptical 
paths each with different radii. Six different SAR system 
configurations were simulated for each terrain feature. 

First, the simple models used to generate the elevation 
data bases characterizing the terrain are presented. Then 
the SAR system parameters are discussed and finally, the ra- 
dar image simulations are presented. 


12.1 SIMULATION OF SIMPLE GEOLOGIC TERRAIN FEATURES 

A prerequisite to performing the system tradeoff study 
using radar image simulation is the modeling of the geologic 
terrain features. These models must then be used to gener- 
ate terrain elevation data for input to the sensor simula- 
tion algorithms. The higher the level of realism in the 
simulated terrain the more useful the results of the sensor 
simulation and ultimately the more worthwhile the tradeoff 
study. However, as a first step only simple terrain fea- 
tures were modeled and simulated. (On-going research is be- 
ing devoted to more realistic simulation of geologic terrain 
features.) Simple geologic features were selected because 
they are easy to implement and to evaluate; also, interpret- 
ing the resulting sensor images is also simple. 

The simple terrain features were constructed by a two 
step process. The first step was to specify the location of 
the peak of the hogback ridges. For the linear ridges this 
was done simply by creating a binary digital image with each 
i,j pixel value assigned 255 if i,j was included as one lo- 
cation of the peak and zero otherwise. For the elliptical 
ridges the i,j pixel values assigned 255 were constrained to 
follow an elliptical path. (See Figures 21 and 22 for the 
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ridge location maps.) The second step in the process was to 
select the cross-section of the ridge which was to follow 
the path specified by the location maps, and to simulate the 
terrain feature by overlaying the selected cross-section or- 
thogonal to the curve defining the peak positions. The 
cross-secti on used for this experiment is shown in Figure 
23. 



Figure 21: Angular Orientation of Linear Hogback Ridges 


Two elevation data bases were thus constructed. However, 
to perform the sensor simulation each terrain location x,y 
(now on the i,j digital matrix position) must have an asso- 
ciated microwave scattering characteristic assigned as well 
as elevation. One microwave scattering characteristic was 
assigned to the entire scene, i.e.* each data base matrix 
location i 5 j has a different elevation but the same scatter- 
ing properties. Figure 24 is the graph of the microwave 
backscatter coefficient versus the angle of incidence used 
here. Note that this ter rain property (scattering) is a 
function of several system parameters, e.g., frequency and 
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Figure 22: Geometry of Elliptical Hogback Ridges 


polarization; thus in any tradeoff analysis the backscatter 
versus angle would change with these system parameters. In 
the simulation algorithms, the local angle of incidence is 
calculated and its corresponding backscatter value is used 
to determine the expected value of the return power for each 
terrain location. 
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Figure 24: Backscatter Coefficient Data 

for Radar Simulation 
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12.2 SIMULATED SAR SYSTEM CONFIGURATIONS 


Six different SAR system configurations were simulated 
for each of the data bases described in the preceding sec- 
tion. Table 3 contains the various system parameters for 
each configuration. This preliminary study was restricted 
to evaluating only two system parameters, radar angle of in- 
cidence and averaging. 


TABLE 3 

System Parameters 



Resolution 

Radar Angle 
of Incidence 

Number of 
Look Averaged 

1 

25m 

23° 

1 

2 

25m 

23° 

4 

3 

25m 

40° 

1 

4 

25m 

o 

O 

4 

5 

25m 

o 

O 

LO 

1 

6 

25m 

o 

O 

<£> 

4 


To interpret the images, the simulation procedure is very 
briefly reviewed. For each data base matrix location i,j 
(in this case the pixel spacing in the data base was set to 
6.25 m) an average returned power is calculated. Incorpo- 
rated into its calculation are all radar and geometric ef- 
fects, e.g., layover, shadow, etc., for the selected angle 
of incidence. A matrix of averaged returned power is thus 
created. This matrix can be scaled and displayed as an im- 
age. The images are totally deterministic in nature and 
have also been generated for each of the three angles of in- 
cidence for each data base. These images allow for the in- 
vestigation of deterministic tradeoff between sensor and 
terrain parameters. Even though these same tradeoffs could 
be determined analyti cal ly, the image analysis approach is 
superior because the ultimate users of the data are gener- 
ally image interpreters. 
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A radar hologram is generated next directly from the 
power matrix by Monte Carlo techniques. Various resolution 
and averaging effects are simulated by applying a linear 
spatially invariant filter. This filter can be implemented 
in either the spatial or spatial frequency domain. The 
equivalent width of the filter's impulse response determines 
the simulated resolution. The output of the filter must be 
formatted for display. The transformation between the out- 
put complex pixel value and the output real pixel value can 
be arbitrarily specified. Common point transformations in- 
clude displaying the pixel power, pixel amplitude on the lo- 
garithm of the pixel power. The SEASAT-A SAR system dis- 
played the pixel amplitude; thus for these simulations the 
pixel amplitude was displayed. Noncoherent integration can 
be achieved either by directly mimicking a discrete mixed 
integration operating in the spatial frequency domain or an 
equivalent continuous mixed integration simulated by a fil- 
ter in the spatial domain. Either approach produces the re- 
quired averaging versus resolution trade-off. 

The images presented in the following section were gener- 
ated by techniques described above. The discrete mixed in- 
tegrator was simulated to produce the four look images. 


12.3 THE SIMULATION RESULTS 


A series of radar simulations were generated for linear 
and elliptical hogback ridges on level ground. The geome- 
tric cross-sections of the features are triangular and sym- 
metrical in general, with an angle of 40 degrees between the 
triangle base and sides. The maximum height of all the fea- 
tures is approximately 50 m as seen in Figure 23. A partic- 
ular radar backscatter cross-section function (which varies 
with angle of incidence) was assigned to both the features 
and the level background (Figure 24). The same statistic 
characteristics were assigned to both the ridges and the 
level ground, i.e., Rayleigh statistics. For each simula- 
tion three images were generated which represent the power 
map of calculated average power based solely upon the radar 
equation, taking into account all considerations 
(e.g., local slope) except fading. The second image of each 
set represents a fully coherent, 25 m by 25 m resolution de- 
piction of the scene. The third image of each set has the 
same resolution, but has four independent samples averaged 
(without sacrifice of resolution). 

In the sensor images of the linear hogback ridges with an 
incidence angle of 23 degrees, there are no shadows in the 
simulation but the large local angles make the backslopes of 
the features appear dark in contrast with the high reflec- 
tivity of the level background (for the no-fading case), as 
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:een in Figure 25. The radar look direction goas from left 
to right. Because of the geomtery of this scene, layover is 
apparent for the first two ridges (counting from the verti- 
cally aligned ridge, left to right). Other functions of ge- 
ometry that are noticeable are the narrowing foreslope rep- 
resentations with increased ridge angle away from vertical 
in the image. The horizontally aligned feature appears to 
have a flat top (and higher reflectivity than the sides of 
the ridge) because of the facet model effects. Figure 25 
(1 look, 25 m resolution) shows the deleterious effects of 
poor si gnal -to-noise ratio characteristic of fully coherent 
sensor images. The horizontal feature has almost disap- 
peared for this noise level. Averaging four looks, as is 
illustrated in the third image of Figure 25, restores the 
detectability of the fifth feature, although its characteri- 
zation in this scene does not hint that the underlying fea- 
ture is a linear ridge. 

When the angle of incidence is increased to 40 degrees 
(altitude unchanged) the background backscatter variation is 
seen upon comparison with Figure 25. Again, local angle ef- 
fects darken the backslopes; there is still no shadow given 
this particular geometry. Significant foreslope compression 
also occurs in this case; foreslope brightening is appar- 
ently due to a zero degree local slope for the vertically 
aligned ridge. The horizontal ridge is somewhat difficult 
to discern even in the "no-noise" case (Figure 26, power 
map). Figure 26 (1 look case) illustrates the loss of de- 
tectability for that ridge which is aligned parallel with 
the radar look direction (horizontally). Because a multi- 
plicative noise model can be used to describe the faded 
scenes, note that the brightness variations are smaller for 
the 1 look case than for the 4 look case in Figure 26. (The 
mean-to-standar d deviation ratios are approximately equal, 
however.) The 4 look case of Figure 25 shows that averaging 
four looks does not really restore the fifth ridge, but does 
improve the i nterpretabi 1 i ty for the four other hogback 
ri dges . 

With an angle of incidence equal to 60 degrees the back- 
ground region has dropped greatly in reflectivity, and sim- 
ple geometric calculations show that shadows now exist for 
^the backslope regions. Again layover can occur and fore- 
'slopes are compressed in this representation; they are very 
bright because of local angle considerations. In the power 
map of Figure 27, the horizontally aligned ridge is all but 
gone, disappears i>~? ‘ihe speckle of the 1 look case, and is 
not visible in the averaged image (4 look case). 

The elliptical hogback ridges were simulated for sensor 
parameters illustrated previously. For the case of 23 de- 
gree angle of incidence. Figure 28, power map case, there is 
no shadow, but very dark (relative to the background) back- 










slopes due to local angle variations. Layover is seen at 
the ground-ridge interface In the image. The bright line 
along the crest of the ridge is due to facet model effects. 
The portion of the r*dge oriented parallel with the radar 
look direction will become less detectable, as in the noisy 
1 look image of Figure 28. Averaging four looks Improves 
the Interpretabl 1 ity of such a ridge geometry (and orienta- 
tion) as viewed in the figure. 

At a 40 degree angle of Incidence, the elliptical ridge 
Is more difficu 1 * to discern (given this particular back- 
scatter model) figure 29 (power map). The darkening is 
again due to loc>_ angle and not to shadow. Portions cf the 
ridge become difficult to see in the faded scene of the 1 
look case of Figure 29, whereas averaging to four looks just 
barely allows the feature to be seen (l.e., the part of It 
which is parallel with the radar look direction). The fact 
that it is a portion of such a ridge type would probably not 
be recognized without supplementary knowledge. 

The 60 degree angle of incidence case is shown (power 
map) in Figure 30. Given this specific combination of back- 
scatter model, ridge model, and orientation, the greater 
part of the geologic feature disappears in the fully coher- 
ent scene of the 1 look case in Figure 30. Averaging four 
looks helps to determine the outline of the feature, but the 
image in Figure 30 provides few clues for the identification 
process . 

Two major conclusions can be stated from these results. 
First, even for a very simple geologic feature, a linear 
hogback ridge, the character of its radar image changes 
drastically with Its orientation to the radar look direc- 
tion. Therefore, it will be difficult if not impossible to 
determine the optimum sensor parameters fo' detecting gen- 
eral linear hogback ridges. However, simulation will pro- 
vide an indication of the geologic mapping potential with 
respect to linear hogback ridges at > specific orientation 
in relation to the radar given different sensor configura- 
tions. Further, these results indicate that totally auto- 
matic detection algorithms for linear hogback features will 
be very difficult to develop because of the diversity of 
their image manifestation. 
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Figure 2 *-. Elliptical Ridges and *0° Angle of Incidence. 
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Figure 30 : 


Elliptical Ridges and 60° Angle of Incidence. 





Chapter 13 
CONCLUSIONS 


The analysis portion of the work reported herein shows 
promise for texture measures of SAR images by the 6LC matrix 
approach despite the poor si gnal -to-noise ratio in the N=4 
SEASAT SAR images employed in this study. The look angle 
dependence of the GLC results must be investigated because 
these effects must be factored into future use of automated 
texture analysis. Point statistics investigations are being 
pursued and will be reported for phase two, along with an 
update on texture analysis. 

The synthesis portion of this work produced geologic fea- 
ture data bases for linear and elliptical hogback ridges. 
The radar image simulations were produced assuming one cate- 
gory (‘'ackscatteri ng type) was present, and clearly demon- 
strated angle of incidence and fading effects. The geologic 
features in this first phase of work are deterministic; ran- 
dom terrain models will be a refinement of the procedure and 
are being pursued for phase two. 
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Chapter 1 
INTRODUCTION 


To date, relatively little research has been devoted to 
the development of optimum specifications for the radar sen- 
sor, or for optimum strategies for information extraction 
from radar data. Careful examination of the problem readily 
explains the situation. Only a limited number of radar sys- 
tems have been available for geologically dedicated research 
programs, and evaluation has consisted of observations de- 
rived from analysis of imagery. Research has been conducted 
with only a limited number of combinations of sensor parame- 
ters such as frequency, polarization, depression angle (or 
angle of incidence), look direction, and resolution. 

This deficiency in control of system parameters is over- 
come with the advent of computer-modeling techniques via ra= 
dar image simulations. Definitive studies to optimize the 
complete radar remote sensing system for geology can now be 
made. Image simulation provides a cost-effective method for 
modeling terrain variation, predicting of results, and de- 
fining of optimum sensor parameters. 

The first phase has recently been completed of a program 
to determine the interpretation strategy and best radar sen- 
sor configuration for geological applications. In this ini- 
tial work, terrain modeling and radar image simulation have 
bee, sed to perform extensive parametric sensitivity stud- 
ies. Several terrain models have been created and numerous 
radar image simulations have been produced from them. It is 
concluded that computer-model i ng via image simulation pro- 
vides a cost-effective method for evaluating terrain varia- 
tions, predicting results, and defining optimum sensor para- 
meters. 
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Chapter 2 
GEOLOGIC MODELS 


In this initial work, several geological terrain models 
have been created. These models represent the textbook case 
in that they are composed of relatively simple geometric 
shapes that are analytically defined, and that they are com- 
pletely deterministic being totally devoid of geologic noise 
( i e . , drainage patterns, weathering effects, surface cover 
boundaries over the landscape model, etc.). They provide, 
however, ideal constructions to evaluate how the radar scene 
changes for a specific landscape model as radar sensor pa- 
rameters are changed. 

An example of such a geological (landscape) model is 
presented in Figure 1. Additional models have been con- 
structed, but there is insufficient space to describe more 
than this one. This model was created mathematically. It 
consists of a grid matrix containing 512 x 512 elements. 
Each element symbolically represents a square region of 
15.625 m on the ground (ie., an area of 244.14 sq m). The 
complete matrix, then, represents a square region of 
64 sq km; each side being 8 km long. Each entry in the ma- 
trix specifies the terrain elevation at a point on the 
ground. The maximum elevation, along the crest-line, was 
set at 100 m. 

As can be seen in the figure, the landscape model illus- 
trated in the three-dimensional plot is that of a breached 
anticline and sync line. The anticline is readily identified 
as the structure having the long nose and infacing escarp- 
ments, while the syncline has no nose and outfacing escarp- 
ments. Slope angles have been set for these features from a 
plunge angle of 8 degrees (on the nose of the anticline) to 
approximately 60 degrees for the very steep escarpments, and 
to approximately 23 degrees on the more gradual slopes. 

This model was created mathematically as a sinusoid ro- 
tated 8 degrees from the surface tangent plane so that it 
plunges into the ground. The period of the sinusoid was set 
so that both the desired slope angles and feature scale 
would be attained. As constructed, the sinusoid produced 
gradual slopes. When the elevation reached 100 m at any lo- 
cation, the sign of the function was reversed and a scale 
factor applied to accelerate the drop back to level ground. 
In this way, gradual slopes rising to a height of 100 m and 
terminating in steep escarpments were created. 
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Chapter 3 

RADAR IMAGE SIMULATION 


A very general approach to simulation has been developed 
at the University of ArkansasCl]. Figure 2 illustrates a 
block diagram of the computer programs. As can be seen from 
the figure, the simulation computer programs can produce se- 
lected intermediate and final results. The required inputs 
to the simulation programs are a data base and a backscatter 
file. Five different intermediate and final results can be 
produced. These are: 1). power map. image; 2). slant-range 
image; 3). ground-range image; 4). slant-range image with 
noise; and 5). ground-range image with noise. Antenna and 
resolution effects can be incorporated via processing either 
in the spatial domain via convolution, or in the spatial 
frequency domain via the Fast-Fourier Transform and multi- 
plication. These are all separately discussed. 
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Figure 2: Simulation Computer Programs 


3.1 DATA BASE 

The data base is a digital replica of the ground and is 
the input mechanism for getting the terrain model into the 
computer simulation programs. This digital representation 
symbolizes the ground and its surface materials (ie., its 
coverl in a grid matrix having at least four dimensions; two 
for the horizontal location of each point, one for its ele- 
vation, and at least one for its cover type. Digital ter- 
rain data (ie., digital representations of terrain eleva- 
tions) are used as the spatial information in a data base. 
These data represent either actual terrain elevations or 
mathematical models of terrain elevations. The data plotted 
in Figure 1 represents one kind of data base. 
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3-2 BACKSCATTER FILE 

The backscatter file contains a model for the reflectiv- 
ity for each terrain cover type present in a data base. It 
is these backscatter data which model the interaction be- 
tween the transmitted electromagnetic energy and each cover 
type. They account for important properties both of the 
ground such as the complex permittivity and surface rough- 
ness, as well as of the radar such as the transmission fre- 
quency and polarization. Empirical backscatter data such as 
those published by UlabyE2], or Cosgriff et a 1 . C 33 are used 
where possible. Theoretical backscatter models such as Fung 
and ChanC43, or HevenorC53 are often important. 


3 . 3 SIMULATION MODEL TRANSFER FUNCTION 

For simulation purposes, a radar is modeled as the system 
illustrated in Figure 3. As can be seen from the figure, 
the transfer function can be symbol! icaly written as 


D * y logCf (h*P r ) ] + K 


( 1 ) 


where D symbolizes the density of silver grains in the film 
of an image, y is the film gamma (ie., the slope of the D vs 
logE curve), and K is a constant depending upon the develop- 
ment processing as discussed by GoodmanC63. This result ar- 
ises from convolving the systems blurring function (h) with 
the input power (Pr). The input power is the average return 
power (Pr) predicted by the radar equation[7,8] multiplied 
by random noise (N 0 ). 
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Figure 3: Transfer Function Model 


3.4 POWER MAP IMAGE 

The power map represents, in imape form, the power den- 
sity being reradiated from each point on the ground back in 
the direction of the radar antema. 'She power map has no 
actual radar analogue, but it is an important interim result 
because it provides diagnostic evidence of backscatter ef- 
fects without the perturbations of added geometric, or pro- 
pagation distortions. 




90 



’ c. * m «** * « a PiJftiM? PC? 

♦ V r - *' 


3.5 SLANT RANGE IMAGE 

The slant-range image can be considered as either an in- 
terim product or a final image. The simulated slant-range 
image is directly analogous to an actual slant-range image. 
It embodies all the geometric distortions inherent in radar. 
The energy being returned to the radar has been sorted and 
re-arranged by time, or distance, in accordance with the ge- 
ometry of the radar and the location of the individual reso- 
lution cell. Geometric distortions such as foreshortening, 
layover and shadow all are incorporated and properly dis- 
played. In essence, the slant-range image represents a map- 
ping from the spatial coordinates (x,y,z) of the ground into 
the image coordinates (r,y); range and azimuth, respec- 
ti vely . 


3.6 GROUND RANGE IMAGE 

The simulated ground-range image is directly analogous to 
an actual radar ground-range image. It represents a nonli- 
near mapping from a slant-range image space back into a 
ground-range (x,y) image space. All the distortions due to 
elevation (relief) in a slant-range image are accentuated 
when mapped into a ground-range image. This mapping pre- 
serves the geometric fidelity, or orthographic arrangement 
of features on a level surface such as the mean ground sur- 
face, but the location of any point above or below the flat 
surface, distorted in the slant-range image, is further dis- 
torted in the ground-range image. 
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Chapter 4 

RESULTS AND DISCUSSIONS 


Simulated radar images have been prepared for several ge- 
ologic models. Simulation is used because it affords pre- 
cise control over both input (ie., the geological model) and 
transformation of the input into a desired visual ouput for- 
mat. Since each step is mathematically controlled, direct 
causal relationships and inferences can be discovered easily 
from the deterministic relationship between input and out- 
put. 

Figures 4-7 present one set of simulated radar images 
created from the geologic landscape model illustrated in 
Figure 1. Figure 4 presents a power map. Figure 5 is a 
slant-range image. Figure 6 is a ground-range image, and 
Figure 7 is identical with Figure 6, except that noise has 
been added. This set of images simulates the results of an 
X-band radar (ie., one operating in the 5-12 Ghz region) 
with the radar being flown to the right across the top of 
each image from an altitude of 800 km. The radar operates 
with HH polarization (ie., Horizontal transmit and Horizon- 
tal receive) and the antenna is symbolically tilted so that 
it illuminates the terrain at an angle of 23 degrees in the 
center of the scene. 

Observe in Figure 4 how well the spatial geometry and de- 
tail are preserved in the power map because none of the pro- 
pagation phenomena inherent in radar have been introduced. 
This means, then, that the image it presents is unrealistic 
because it is unattainable from an actual radar. It con- 
tains a significant amount of information for our purposes, 
however, because it shows the normalized power density being 
reradiated back toward the antenna. Loosely stated, it is 
the photographic alternative to a radar image; as if it were 
somehow possible to illuminate the ground at X-band and pho- 
tograph it directly. 

The image presented in Figure 5 looks very different from 
the one in Figure 4. This is because the image in Figure 5 
is a simulated slant-range radar image whereas the image in 
Figure 4 is not even a radar image. The slant-range image 
is a natural consequence of ucing radar to produce images. 
A ground distance of 8 km out to the side of the radar has 
been mapped (compressed), at the proper scale, into a 
slant-range distance of 3.13 km. We see at this angle an 
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Figure 12: 60 Degree 
Power Map Image 


Figure 8: 40 Degree 
Power Map Image 


Figure 4: 23 Degree 
Power Map Image 


Figure 13: 60 Degree 
Slant Range Image 


Figure 9: 40 Degree 
Slant Range Image 


Figure 3: 23 Degree 
Slant Range Image 


Figure 1^ : 60 Degree 
Ground Range Image 


Figure 10: 40 Degree 
Ground Range Im*ge 


Figure 6: 23 Degree 
Ground Range Image 


Figure 15: 60 Degree 
Ground Range Image 
with Noise 


Figure 11: 40 Degree 
Ground Range Image 
with Noise 


Figure 7: 23 Degree 
Ground Range Image 
with Noise 
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extreme example of compression, foreshortening, and layover. 
Since the data base is mathematical and, thus, smooth, it is 
easy for us to observe these effects. 

Figure 6 shows a ground-range version of the slant-range 
image. In the ground-range image the compression evident in 
Figure 5 has been removed. Note how spatial distortion is 
exaggerated for features above (or equally for features be- 
low) level ground. 

Figure 7 illustrates the ground-range image with statis- 
tical noise incorporated* Even for this simple, smooth, 
well-behaved figure, it is surprising to see how noise makes 
interpretation more difficult and less certain. If we had a 
feature contaminated with geologic noi se (ie., small scale 
perturbations such as drainage and weathering effects, vari- 
ations in the rock types, etc.), it is easy to imagine how 
difficult it would be to assess all of these effects, sepa- 
rately. 

Figures 8-11 present the same sequence of simulated im- 
ages from the same data base, but for a different angle. 
The angle for this set is 40 degrees. Evaluation of this 
set of images, like the previous one, shows that layover is 
still present though less pronounced. Spatial distortion, 
though, is much less severe. Note, however, how effectively 
the shape information is masked in the presence of noise. 

Figures 12-15 present another sequence of simulated radar 
images for an angle of 60 degrees. Note the introduction 
of shadow. In this set of images, propagation effects are 
seen to be minimal and it is readily observed that the 
ground-range image agrees very favorably with the power map. 
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Chapter 5 
CONCLUSIONS 


These preliminary results show that digital image simula- 
tion and computer-modeling afford important advances for 
this kind of study. They provide cost-effective methods for 
evaluating terrain variations and sensor parameter changes, 
for predicting results, and for defining the optimum sensor 
parameters for a geologic application. In addition, they 
provide extremely useful mechanisms for teaching radar geol- 
ogy, and for defining an optimum interpretation strategy for 
a geologic application. 
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Chapter 1 
INTRODUCTION 


Variation of the intensity or shade of grey (ie., tone) 
in radar images portrays the relative strength of the radar 
signal returned to the receiver from point to point on the 
ground. Variation of the radar signal arises from numerous 
interacting and complex causes. Two fundamental causes ar- 
ise; first the interaction between the ground and the trans- 
mitted electromagnetic energy illuminating the ground, and 
second the geometric, or propagation, phenomena associated 
with the fact that radar is a ranging device (ie., it re- 
cords the returned signal versus time and, thus, orders fea- 
tures according to their individual distances, or ranges to 
the radar). 

Typically, the percentage of energy backscattered from 
slopes facing a radar is greater than from level terrain and 
from slopes facing away from a radar. Foreslopes are gener- 
ally depicted in an image as brighter tones than level 
ground, and brighter still than backslopes. For the various 
terrain cover types present, this is governed by the back- 
scatter trends versus the local angle of incidence (ie., the 
angle between the local vertical and a line to the radar). 
The geometric relationships are illustrated in Figure 1. 

Geometric distortions arising from propagation phenomena 
perturb this picture. Typically in an image, areas of for- 
eshortening and layover are represented by bright tones and 
regions of shadow by dark tones. This is because areas suf- 
fering foreshortening and layover cause the return signal 
from a relatively large ground area to be mapped into a 
smaller one. At the same time, shadowed regions are areas 
for which very little, if any, signal is returned to the ra- 
d ar . 

In a radar image backscatter and propagation effects are 
simultaneously present and impossible to separate. Only 
through the experience gained in using radar imagery is it 
feasible to estimate the relative contribution of each ef- 
fect for a specific tone in an image. Where it has been 
possible to check, we have observed a tendency for users to 
overestimate the backscatter contribution. 

One unique method for observing the relative contrib- 
utions of tkese effects is afforded by radar image simula- 
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Figure 1: Radar Geometry 


tlon. Digital terrain data are used to model radar Image 
formation. Both backscatter and propagation effects are 
separately modeled. These are serially incorporated and the 
Image expression of each Is noted. Sequences of images are 
presented Illustrating these effects across a range of 
slopes and angles of incidence. The conclusions reached are 
that at angles of incidence smaller than the average slope 
of the terrain in a region, propagation phenomena dominate. 
As the angle of incidence Increases beyond this, the radar 
image portrays an increasingly more faithful representation 
of the backscatter from the ground. It is also show” that 
digital simulation affords an important tool for evaluating 
complex interactions between the ground and radar, for 
training users in radar image Interpretation, or for select- 
ing optimum sensor parameters for specific applications. 



Chapter 2 

PROPAGATION PHENOMENA 


An old radar image (circa 1957) of the Aichilik River 
Area in Alaska is presented in Figure 2 to illustrate the 
visual impact created by the various propagation phenomena. 
This radar image is a ground-range presentation recorded by 
the AN/APQ-56 system. 

Figure 3 presents a topographic map for the region por- 
trayed in the radar image. The topographic map illustrates 
that the mountainous terrain is relatively homogeneous 
whereas the radar image does not give that impression at 
all. The average slope of the terrain in the area is ap- 
proximately 30 degrees. As can be seen in the image, the 
line A-C-E splits the image between the region in layover 
and the region not in layover. The line corresponds to an 
angle of incidence of approximately 30 degrees. Features 
above this line and, thus, nearer the radar, are in layover, 
while those below the line are not. 

In the middle of the image (ie., in the mid-range region) 
note that foreslopes are brighter than backslopes, as ex- 
pected. Foreshortening is maximum for features near the 
line A-C-E with a complete foreslope being mapped into a 
single line in the image. As the angle of incidence is in- 
creased and we move down the image away from the line A-C-E, 
note that foreshortening decreases and foreslopes are pro- 
gressively mapped into larger areas in the image (closer to 
the actual scale size). In the far-range region, shadow be- 
gins to dominate the image presentation, foreshortening is 
minimized, and image tones are very nearly representative of 
the percentage of energy backscattered to the radar. 

At the top of the image (ie., in the near-range region) 
layover dominates for mountainous terrain and it is very 
difficult to interpret the image. The VEE's so vivid in 
this region are the image expression of foreslopes in lay- 
over being exaggerated by a ground-range display. The tip 
of a VEE represents the topmost part of a mountain and is, 
thus, imaged closest to the radar. 
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Chapter 3 
BACKSCATTER 


Radar return from the ground is described by a (the dif- 
ferential scattering cross-section, or scattering coeffi- 
cient)!^]. The scattering coefficient, a 0 , determines the 
electromagnetic energy reradiated back to the radar from the 
resolution cell. A resolution cell is defined for a short- 
pulse, side-looking, imaging radar to be the instantaneous 
area on the ground from which backscattered power arrives 
simultaneously at the radar antenna. A resolution cell is 
considered to consist of many individual features or objects 
relative to the scale of the radar wavelength, and each of 
these is termed a scattering element. 

Figure 4 illustrates the trends of o° versus local angle 
of incidence (o p) for 4 typical surfaces: 1). specular or 
smooth; 2). slightly rough; 3). moderately rough; and 4). 
extremely rough. The scale of the roughness in terms of the 
wavelength of the incident radiation is what determines this 
characteristic. A given surface may appear almost smooth at 
sufficiently low frequencies (ie., at long wavelengths) and, 
thus, exhibit nearly specular reflection. As the frequency 
is increased (ie., as wavelengths become shorter) the same 
surface will appear progressively rougher until it becomes 
extremely rough and, thus, exhibits nearly uniform backscat- 
ter£23. The scale of roughness of a surface relative to the 
radar frequency, within the limitations of this simplistic, 
general discussion, defines the shape of the vs curve. 
As can be seen, the surface that is smooth relative to the 
radar wavelength produces an extremely strung return only 
when the local angle of incidence is zero degrees with the 
strength of the return falling off rapidly as the angle is 
increased. An extremely rough surface relative to the radar 
wavelength is shown to have a ct° vs 0£ trend which is almost 
a horizontal line. In addition, a 0 vs 0£ trends are shown 
for two intermediate surfaces. As might be expected, the 
slightly rough surface exhibits a more specular trend 
whereas the moderately rough one exhibits a more uniform 
trend. 

Generally speaking, we can summarize backscatter trends 
In imagery as follows. For a given scale of roughness and 
complex permittivity, the smaller the local angle of inci- 
dence the higher the backscatter and the brighter the image 
tone. Foreslopes, typically, have small angles (o^< 30 de- 
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Figure 4: Backscatter Trends 


grees) and, thus, relatively high backscatter. Couple this 
with the geometric effects of compression and layover and 
the image tone is very bright; often it is saturated 
(ie., white). Once recorded in imagery, the two effects, 
backscatter and geometric, cannot be separated. In addi- 
tion, the resulting image is probably uni nterpretabl e in 
such areas due to geometric distortions. The greater the 
0 £, as is frequently the case for backslopes, the darker the 
image tone. Within 5 degrees to 10 degrees of grazing 
(ie., 80 degrees <0£<9O degrees) the backscatter begins to 
fall off very rapidly. 
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Chapter 4 

RESULTS AND DISCUSSIONS 


Computer-generated radar images have been prepared via 
the University of Arkansas' radar image simulation pro- 
grams£3] to illustrate the separate effects caused by propa- 
gation phenomena and backscatter. Five different intermedi- 
ate and final results can be produced. These are: 
1). power map; 2). slant-range image; 3). ground-range im- 
age; 4). slant-range image with noise; and 5). ground-range 
image with noise. Simulation is especially beneficial for 
separating propagation effects from those of backscatter be- 
cause these effects, inseparable in actual imagery, can be 
modeled as serial events and the impact of each can be as- 
sessed, separately. 

Figures 5-8 present one set of simulated radar images 
created from an input which is deterministic and ideally 
suited to such a study. It is the geological landscape 
model of a breached anticline and syncline structure previ- 
ously presented[33 . This particular geologic model was cre- 
ated mathematically. As can be seen in Figure 5, the model 
represents a textbook example of an anticline and a syncline 
structure. The anti cline is readily identified as the 
structure having the long nose and infacing escarpements , 
while the syncline has no nose and outfacing escarpements. 

Figure 5 pre-ents a power map. Figure 6 is a slant-range 
image. Figure 7 is a ground-range image, and Figure 8 is 
identical to Figure 7, except that noise has been added. 
This set of images simulates the results of an X-band radar 
(ie., one operating in the 5-12 GHz region) with the radar 
being flown to the right across the top of each image at an 
altitude of 800 km. The radar operates with HH polarization 
(ie., Horizontal transmit and Horizontal receive) and the 
antenna is symbolically tilted so that it illuminates the 
terrain at an angle of 23 degrees in the middle of the 
scene . 

Observe in Figure 5 how well the spatial geometry and de- 
tail are preserved in the power map. This is because in the 
power map none of the propagation phenomena inherent in ra- 
dar have been introduced. This means that the image it 
presents is unrealistic because it is unattainable from an 
actual radar. It contains a significant amount of informa- 
tion for our purposes, however, because it shows the normal- 
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Figure 5: 23 Degree Figure 9: 40 Degree Figure 11: 70 Degree 

Power Map Image Power Map Image Power Map Image 



Figure 6: 23 Degree Figure 10: 40 Degree Figure 14: 70 Degree 

Slant Range Image Slant Range Image Slant Range Image 
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Figure 8: 23 Degree 
Ground Range Image 
with Noise 


Figure 12: 40 Degree 
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ized power density (ie., it is proportional to the actual 
backscatter) being reradiated back toward the antenna. 
Loosely stated, it is the photographic alternative to a ra- 
dar image; as if it were somehow possible to illuminate the 
ground at X-band and photograph it directly. 

Note in Figure 5 how each slope facing the radar (fore- 
slope) is bright except the very steep ones, and each slope 
facing away (backslope) is dark relative to the tone of the 
background. Thus, as would be expected from typical back- 
scatter trends (eg., Figure 4), foreslopes are brighter and 
backslopes are darker than the background. 

The image presented in Figure 6 looks very different from 
the one in Figure 5. This is because the image in Figure 6 
is a simulated slant-range radar image whereas the image in 
Figure 5 is not even a radar image. The slant-range image 
is a natural consequence of using radar to produce images. 
A ground distance out to the side of the radar of 8 km has 
been mapped (compressed), at the proper scale, into a 
slant-range distance of 3.13 km. We see at this angle an 
extreme example of compression, foreshortening, and layover. 
Follow the trace of the base of the geologic structure 
across the image. It is easily identified throughout most 
of the image as the boundary between the background and a 
darker tone. This is especially visible for the ridges run- 
ning across the central portion of the image swath 
(ie., from the leftmost curve to the rightmost one). After 
finding the trace of the base, compare with Figure 5 and ob- 
serve how the top of the ridge has been displaced toward the 
radar (ie., toward the top of the image). This is easily 
ascertained by finding the trace of the top of the ridge as 
the boundary between a brighter tone and the background. 

Figure 7 shows a ground-range version of the slant-range 
image. In the ground-range image the compression evident in 
Figure 6 has been removed, and the trace of the base of the 
geologic structure preserves spatial geometry. Note how 
spatial distortions are exaggerated for features above (or 
equally for features below) level ground. 

Find, now, the trace of the top of the ridge in this fig- 
ure (7). It is easy in this figure to trace the top across 
the image. Observe how in the central portion of the swath 
the top is imaged before the bottom and is, thus, laid 
over . As the trace of the top of the ridge is followed 
around the curve going down the image, two crossings of the 
base occur. This represents a transition in the local angle 
of incidence from a region of layover to one of foreshorten- 
ing and back to layover again. Following the trace of the 
top of the ridge up the image shows similar effects. 
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Again, compare the ground-range and power map images 
(Figures 5 and 7). Observe the spatial distortions that are 
evident, Observe also how tones are distorted. Slopes that 
are in layover are significantly brighter than predicted in 
the power map image. It is interesting to evaluate only the 
central ridge. The true spatial geometry is depicted in 
Figure 5. Compare this to the backslope in Figure 7. Since 
the terrain is continuous, the backslope must be the area 
occupying the region in the image from the top of the ridge 
to the base on the bottom side. We know the top of the 
ridge is in layover in this region ( i e . , top imaged before 
the bottom), so the backslope must extend from the top of 
the white trace (ie., top of the ridge) to the bottom of the 
dark trace (ie., bottom on backslope side). The backslope 
is, thus, stretched out and banded. The first band, the 
lower one, is the correct tone for the backslope, as can be 
verified by comparing with Figure 5. The second band, the 
upper bright one, consists of the power returned from part 
of the backslope, the foreslope, and level terrain in front 
of the ridge. In fact, it is only because we have such a 
well-behaved feature that it's fairly easy to interpret it 
properly. 

Figure 8 illustrates the ground-range image with statis- 
tical noise incorporated. Even for this simple, smooth, 
well-behaved figure, it is surprising to see how noise makes 
the interpretation of this image more difficult and less 
certain. If we had a feature contaminated with geologic 
noise (ie., small scale perturbations such as drainage and 
weathering effects, variations in the rock types, etc.), it 
is easy to imagine how difficult it would be to assess all 
of these effects. We see, therefore, how radar image simu- 
lation is valuable for evaluating the separate impacts of 
various phenomena. The basic conclusion drawn from this set 
of images representing a small incidence angle is that pro- 
pagation effects dominate backscatter ones. 

Figures 9-12 present the same sequence of simulated im- 
ages from the same data base, but for a 40 degree angle of 
incidence; power map, slant-range, ground-range, and 
ground-range with noise images, respectively. As can be 
seen from the various images, compression is less severe 
(ie., 8 km mapped into 5.142 km) as are all the propagation 
effects. Evaluation of this set of images, like the previ- 
ous set, shows that layover is 'still present though less 
pronounced. Thus, propagation phenomena dominate backscat- 
ter effects even at this intermediate angle. Spatial dis- 
tortion, though, is much less severe. Note, however, how 
effectively the shape information is masked in the presence 
of noise. 

Figures 13-16 present another sequence for an angle of 
incidence of 70 degrees. Note the introduction of shadow. 
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This can be seen both in the images and by the fact that the 
local angle of incidence has exceeded the grazing angle. In 
the set of images, propagation effects are seen to be mini- 
mal (1e., 8 km mapped into 7.52 km) and it is readily ob- 
served that the ground-range image agrees very favorably 
with the power map. For this large angle set, then, it is 
readily seen that spatial geometry is nearly preserved and 
the image almost records the true backscattered signal rela- 
ti onships. 
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Chapter 5 
CONCLUSIONS 


Simulation was used to decouple image effects thereby al- 
lowing evaluation of the image expression of each separate 
effect. It was shown that propagation phenomena are impor- 
tant distortions, at least for radar imagery collected at 
angles of incidence less than 40 degrees, which must be 
properly accounted for when interpreting such imagery. Im- 
proper allowance for propagation distortions will lead to 
incorrect interpretations and erroneous beliefs regarding 
the properties of the ground because the relative amount of 
backscatter and the terrain slope will likely be misjudged. 
From the evidence presented, it was shown that the dominat- 
ing factor in determining the level of brightness in an im- 
age is propagation phenomena for small angles of incidence, 
and backscatter for large angles. From this study it is ob- 
vious that digital simulation affords an important tool for 
correctly interpreting image distortions, and for training 
users in radar image interpretation. 
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Chapter 1 
INTRODUCTION 


It has been held, traditionally, that the value of radar 
imagery is increased for most applications when the angle of 
incidence is large. Figure 1 illustrates the geometry for 
an imaging radar over a curved earth and identifies various, 
important relationships such as angle of incidence. With 
the advent of spaceborne imaging radars such as the SEASAT 
synthetic aperture radar, it is often infeasible to attain 
these desired angles for practical reasons. The tradeoffs 
required between frequency, antenna size, and power availa- 
bility are the most limiting of these reasons. 

Because the large angles which are desired (ie., e>60 de- 
grees) for most geological work are unattainable, a study is 
being conducted to assess the usefulness of imagery col- 
lected at smaller angles for satisfying geological needs. 
The study consists of obtaining and evaluating radar imagery 
collected at numerous angles of incidence for specific geo- 
logical features. Since it is generally impossible to fhd 
suitable imagery of a specific feature illuminated across 
the complete range of angles, radar image simulation is be- 
ing employed to accomplish this. The value of using simula- 
tion and computer-generated imagery is clearly established 
in studies such as this one. 

Preliminary results ini cate that the usefulness of radar 
imagery is highly dependent upon the application. Imagery 
collected at small angles of incidence was found to be use- 
ful for applications involving study of low relief terrain 
because typical trends of the backscatter coefficient versus 
the 1 ocal angle of incidence results in enhanced contrast 
for such features in imagery; a conclusion also reached by 
Ford Cl].. Figure 2 illustrates a curve showing typical back- 
scatter trends. These typical trends are presented for 
their relative shape information, so neither amplitude scale 
nor frequency and polarization are identified. 

Imagery collected at large angles of incidence was found 
to be most useful for applications involving study of moun- 
tainous terrain because of the out! i ni ng effect of radar 
shadow and because geometric distortions are minimized. Im- 
agery collected at intermediate angles was found to serve 
neither type of application best, but probably offers a com- 
promise for viewing mixed terrain types since large angles 
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Figure 1: Radar Geometry 


are difficult to attain. A system having at least two 
selectable angles of incidence, one small and one large, 
would probably be the best one for imaging mixed terrain. 
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Figure 2: tfackscatter Trends 
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Chapter 2 


INCIDENCE ANGLE EFFECTS IN RADAR IMAGERY 


We have observed that the detail and fidelity portrayed 
in radar images of different kinds of terrain is dependent 
upon both the relative relief of the terrain as well as the 
radar angle of incidence. This observation was drawn from 
analyzing computer-generated radar images portraying a sin- 
gle scene imaged at numerous different angles of incidence, 
and it was validated by comparing these to actual imagery. 

The computer-generated images were created by using the 
University of Arkansas' radar image simulation computer pro- 
grams'^]. Four of these simulated radar images are pres- 
ented as ground-range images in Figures 3-6. Each of the 
figures illustrates the same scene imaged from spacecraft 
altitudes at a specific radar angle of incidence; 23 de- 
grees, 40 degrees, 60 degrees, and 80 degrees, respectively. 
The relative relief for every point in the scene was ob- 
tained from a digital terrain tape of the NJ 16-12E standard 
USGS topographic map[4]. Data on this tape recorded the el- 
evation of the area on a 91.44 m sampling grid. 

Each image consists of an array of 512 x 512 points, or 
pixels (ie., picture elements). Each pixel represents a 
ground spot of 91.44 m x 91.44 m. The scene in each figure 
portrays, therefore, a ground-range image of 2192 sq km 
(ie., 46.8 km on each side) in the southern Appalachian ba- 
sin region of Tennessee. Two types of terrain are present 
in this region: 1). maturely dissected mountains, and 
2). elongated mountain slopes. 

Figure 3 presents a computer-generated ground-range image 
of the Tennessee terrain for a spaceborne radar orbiting 
the earth at 800 km and viewing the scene at an incidence 
angle of 23 degrees. In this simulated scenario, the angle 
of incidence is 21.6 degrees at the near-range edge, 23.0 
degrees in the middle, and 24.4 degrees at the far-range 
edge of the scene. Observe in this scene how the maturely 
dissected mountainous terrain is distorted by extreme lay- 
over and further distorted by the ground-range presentation 
resulting in the characteristic VEE 1 s pointing toward the 
radar, as previously di scussed[57T The presentation of 
mountainous terrain in this angular range severely distorts 
the geometric fidelity, thereby making such radar images 
only marginally useful for many geological investigations. 
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Figure 3s 23 Degree 
Ground Range Image 
with Noise 


Figure 4: 40 Degree 
Ground Range Image 
with Noise 



Figure 6: 80 Degree 
Ground Range Image 
with Noise 



Figure 3: 60 Degree 
Ground Range Image 
with Noise 
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Note the sharp detail (ie,, enhanced contrast) and 
apparent fidelity of the presentation in this image of the 
low relief, elongated ridge terrain. Even though compres- 
sion and layover are present, the images of this terrain- 
type is useful for many investigations. The sharpness in 
this area is attributed to the sensitivity of the typial 
backscatter coefficient to small angular changes near nadir. 

A curve showing the typical trends of backscatter coeffi- 
cient versus local angle of incidence is presented in Fig- 
ure 2. Since the backscatter coefficient determines the 
amount of electromagnetic energy reradiated from each reso- 
lution cell back to the radar, it is the fundamental mecha- 
nism which determines the shade of grey (or tone) of each 
pixel. As can be seen in the figure, small angular changes 
around an average angle of 23 degrees result in relatively 
large changes in the backscatter coefficient. The large 
backscatter changes are recorded in imagery as large varia- 
tions in the tone from pixel to pixel. In this way, the ra- 
dar expression of terrain that is very rough but not moun- 
tainous relative to the radar wavelength is recorded with 
enhanced contrast in imagery collected at an angle of inci- 
dence of 30 degrees (approximately), or less. 

In the 30 degrees to 70 degrees region (approximately),, 
the typical backscatter trend shows a marked reduction in 
the si ope . We would expect small angular variations associ- 
ated with very rough terrain when viewed in this region to 
result in imagery having very little contrast. This is pre- 
cisely what is shown in Figures 4 and 5. In Figure 4, the 
angle of incidence is 39.0 degrees at the near-range edge, 

40.0 degrees in the middle, and 41.0 degrees at the far- 
range edge. In Figure 5, these angles are 59.6 degrees, 

60.0 degrees, 60.4 degrees, respectively. Note in both fig- 
ures how the contrast has been reduced for the very rough 
terrain to the point where, in comparison to the presenta- 
tion in Figure 3, these images are no longer indicative of 
the topography. 

The mountainous terrain exhibits only a slight amount of 
layover in Figure 4, but compression is still se- 
vere(ie., 36%). This causes the linear extent of foreslopes 
to be foreshortened, or compressed in an image. The image 
of such foreslopes is largely uninterpretable. Figure 4, 
thus, illustrates a radar image that is only marginally use- 
ful . 

Very little, if any, of the mountainous terrain is laid- 
over in Figure 5. Even though compression still exists at a 
60 degree angle of incidence (ie., 13%), it does not cause 
significant distortions for most terrain. For this terrain, 
the appearance of radar shadow aids interpretation in the 
mountainous region. This image, then, marks the transition 
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point (approximately) from which radar imagery of 
mountainous regions is most interpretable. Despite the ap- 
pearance of shadow, the image of the low relief terrain is 
not nearly as useful as the one in Figure 3. 

The simulated image presented in Figure 6 was produced to 
represent a 79.9 degree angle of incidence at the near- 
range, 80.0 degrees in the middle, and 80.1 degrees at the 
far-range edge. This is a spectacular image even though 
dominated by shadow. It is easy to see why imagery near 
this incidence angle is so highly desired by geologists. It 
is shown only for comparison since a basic premise underly- 
ing our work here is that large angles are presently unat- 
tainable from spacecraft altitudes. Further analysis of 
this image indicates this angular region (ie., near grazing) 
marks (approximately) the point where a typical backscatter 
curve again enjoys a significant slope ( eg, , Figure 2). 
This is a region, just like the one near nadir, where small 
angular changes result in increased image contrast. For 
certain types of terrain, this occurs. Often, though, 
shadow dominates, causing loss of much information, but ra- 
dar imagery in this angular region seems highly desirable 
for the great clarity of presentation it affords for most 
terrain types. 
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Chapter 3 
CONCLUSIONS 


The detail and fidelity portrayed in radar imagery of 
different kinds of terrain is dependent upon both the rela- 
tive relief of the terrain as well as the radar angle of in- 
cidence. The usefulness of radar imagery for geological ap- 
pl i cati ons is , thus, dependent upon the application. 
Analysis of our preliminary results show, even though only a 
very limited range of parameters have been evaluated, that 
large incidence angle radar imagery is best for many geo- 
logic applications. If we accept the premise that these an- 
gles are unattainable, then our preliminary results indicate 
that applications which involve study of low relief terrain 
can be served best by imagery collected at small angles-of- 
incidence (ie., e<30 degrees), and applications which re- 
quire study of mountainous terrain can be served best by im- 
agery produced at the largest angl es-of-i ncidence attainable 
(ie.,e>60 degrees). The angular range in between these 
serves neither low relief nor mountainous terrain as well, 
but offers a compromise. Further study is needed to evalu- 
ate additional parameters to see if a particular combination 
will provide imagery having increased contrast and clarity 
in this intermediate angular range. In general from our 
preliminary work, it would appear that for a mix of terrain 
types the best compromise is offered by a system having at 
least two selectable angl es-of-i ncidence; one small and one 
large. The next best compromise is probably offered by a 
system viewing the scene only at the largest angle-of-inci- 
dence which can be attained within the spacecraft limita- 
tions. 
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